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1.0 SUMMARY 


This program has investigated four-wave mixing (4-wm) in nonlinear waveguides, 
with particular emphasis on waveguides with gain and a resonant enhancement. 

In 4-wm a signal or probe beam is scattered off of a grating induced in a 
nonlinear medium due to interaction with two pump beams. Under appropriate 
conditions, the scattered beam will be an amplified phase-conjugate (PC) or 
frequency-shifted replica of the signal beam. Such a scattered beam has 
applications in laser device phase locking, beam steering, and optical comput¬ 
ing. Phase conjugation by 4-wm in nonlinear waveguide with gain and/or a 
resonant enhancement has the potential to be a high-speed process requiring 
only modest power levels. 

To understand the basic physical processes which determine the formation, 
duration, tunability, and fidelity of the scattered optical signal, a general 
analytical framework has been developed for studying nonlinear optical 
processes in waveguides based on the coupled-mode formalism. This approach 
has the advantage of developing a model of 4-wm which is similar to the 
coupled-wave formalism used to describe nonlinear optics in bulk media. The 
influence of diffusion on those media where the nonlinearity is dominated by 
mobile particles, such as free carriers in a semiconductor under current 
injection, can be incorporated through a Lorentzian dependence on spatial and 
temporal frequency of the interacting modes. Our theoretical analysis shows 
that there are strict limits on the speed and spatial fidelity of the phase 
conjugation process in a multimode waveguide. However, these limits can be 
relaxed by using an array of single-mode waveguides. 

The theoretical framework gives as a requirement for maximum possible fidelity 
in conjugate generation, that there be no degeneracy among the guided modes. 

An experiment using an optical fiber showed that the degeneracy between the 
two orthogonal polarizations in the circularly symmetric medium leads to 
polarization scrambling through both linear inhomogeneities as well as non¬ 
linear interaction. Polarization scrambling reduces fidelity and complicates 
the analysis of the nonlinear interaction. 
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Laser diodes are a particularly interesting nonlinear waveguide due to the 
simultaneous combination of strong nonlinearity, resonant gain, and a Fabry- 
Perot cavity. Previous work has shown conjugate signal generation with 
intensities 1,000-fold larger than the input signal intensity. It has been 
demonstrated that the action of a nearly-degenerate, weak signal beam is to 
phase modulate the output of a laser diode through 4-wm. To be a phase 
modulation, the 4-wm signal must be essentially equal in amplitude with the 
amplified input signal and maintain a specific phase relationship with respect 
to the input signal and the counterpropagating pump beams. In the semicon¬ 
ductor medium under current injection, the nonlinearity acts to minimize 
carrier fluctuations induced by the probe beam through the generation of a 
conjugate beam, oppositely shifted in frequency with respect to the 
free-running diode. 
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2.0 INTRODUCTION 


Nonlinear optical processes using 4-wm in bulk media have been demonstrated 
and studied in a variety of applications (Ref. 1). P.ecently, much of the 
scientific work has gone into the investigation of photorefractive materials 
such as BaTi0 3 . Four-wave mixing in photorefractive materials can be achieved 
at low intensities relative to media using other nonresonant properties or 
stimulated Raman or Brillouin scattering (Ref. 1). There have been a variety 
of studies showing application of 4-wm in BaTi0 3 to the multiple coupling and 
locking of lasers (Refs. 2, 3, and 4), beam steering (Ref. 5), and optical 
computing (Refs. 6 and 7). However, there are stringent bounds on the speed 
of photorefractive materials at low laser intensity (Ref. 8) due to the time 
required for the photo-generated carriers to achieve sufficient density and 
separation to create the internal electric field. Since many applications 
demand high speed and modest power levels, there are limitations on the 
utility of bulk materials in 4-wm applications. 

Phase conjugation by 4-wm in waveguides containing a nonlinear medium with 
gain and/or a resonantly enhanced nonlinearity would not be limited by these 
considerations. Four-wave mixing in waveguides was originally investigated to 
take advantage of the waveguide’s capability to maintain high intensity over a 
long interaction length and the guide's reduced alignment requirements 
(Refs. 9-12). Following these initial reports, there was no further activity 
for several years because there were no suitable nonlinear materials 
(Ref. 13). Recently, there has been a report of 4-wm in a semiconductor laser 
(Ref. 14) and an observation of large optical nonlinearities in GaAs/AiGaAs 
multiple quantum wells (MQW) (Ref. 15). Subsequent analysis of the 4-wm in 
the semiconductor laser has not included the effects of the guided-wave 
structure (Refs. 16-19) and there have been no studies of 4-wm using the MQW 
excitonic resonance in waveguides. A theoretical study is proposed to make 
the connection between the new semiconductor materials, with gain and resonant 
enhancement, and the particular requirements and advantages imposed by multi- 
mode waveguides, as outlined by Hellwarth (Ref. 11), on the problem of 4-wm. 
The goal of this study is to understand the basic physics involved, with 
attention to the feasibility of the system to applications such as laser 
coupling and phase locking. 
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2.1 BACKGROUND 


Hellwarth's article was the first to show the general utility of "generation 
of a PC of a free wave at entrance of a waveguide" (Ref. 11). The system 
is depicted schematically in Figure 1. His analysis pointed out that as long 
as the pump and probe beams occupy many modes and are well guided, the output 
of the 4-wm process will be the PC beam. Further, beam alignment is not 
critical; it is only necessary to inject the beams into the guiding modes. 
Simultaneously, experimental work demonstrated the advantages of the long path 
length, high intensity interaction region possible in waveguides (Refs. 10 
and 12). Lower pump levels relative to bulk media were possible, but the 
ratio of conjugate beam power to input beam power was less than unity. The 
waveguide medium, CS 2 , did not have•a sufficiently strong nonlinearity to pro¬ 
duce stronger outputs. For this reason, the theoretical analysis was limited 
to the regime of small conjugate reflectivities. Subsequent work showed 
equally small conjugate reflectivities in waveguides where CS 2 was used as the 
cladding, rather than guiding, medium (Ref. 20). 

The recent observation of 4-wm in an AiGaAs semiconductor laser showed 
conjugate reflectivities approaching 10* due to the resonant enhancement and 
gain in the laser medium (Ref. 14). The nonlinearity, created by carrier- 
density modulated index gratings, has a response time governed by carrier 
spontaneous lifetimes, independent of excitation intensity (Ref. 18). These 
lifetimes are on the order of 1 ns in III-V semiconductors. Recent calcula¬ 
tions show that response times as short as 10 ps are possible in semiconductor 
lasers (Ref. 19). The theoretical modeling of the single-mode laser diode 
structures has focused on the nonlinear semiconductor medium without consider¬ 
ing the waveguiding aspects of the interaction (Refs. 16-19). Therefore, the 
connection remained to be made between these recent materials studies and the 
issues of 4-wm in multimode waveguides raised by Hellwarth. 

This program has addressed the issue of conjugate generation through 4-wm in 
waveguides containing both nonlinearity and resonant amplification. A general 
formalism for studying waveguide nonlinear optics is developed in Section 3.0. 
The effects of carrier diffusion are analyzed in Section 4.0. These topics 
have been summarized in two publications (Refs. 21 and 22). Section 5.0 
reviews some of the factors which limit the speed and fidelity of conjugate 
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Figure 1. Schematic illustrations of two arrangements for generating an optical 
beam, E, which is PC to an input beam by 4-wm with pump beams 
G and H (Ref. 11). a) Three-beam configuration with F being phase 
conjugated. If the waveguide is a laser, G and H are the intracavity 
oscillating beams (Ref. 14). b) Two beam configuration where one 
of the pump beams, G, is phase conjugated. 


formation in a nonlinear waveguide. An array of multimode waveguides is 
proposed as a way to overcome these limits (Ref. 23). The array format shows 
that the nonlinearity in laser diodes can be used for conjugate formation. 

The results of experimental investigations of 4-wra in laser diodes are 
described in Section 6.0, along with experiments on optical fibers. The laser 
diode work shows that the action of an injected, nearly-degenerate signal is 
to phase modulate the optical output in a way that works to minimize the fluc¬ 
tuations of free carriers in the device. The fiber work illustrates some 
complicating effects of polarization on the nonlinear interaction in a system 
where the two orthogonal polarizations yield nearly degenerate modes. 

Section 7.0 summarizes the results of the program and suggests directions for 
future research. 
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3.0 COUPLED-MODE THEORY FOR NONLINEAR WAVEGUIDES 


3.1 INTRODUCTION 

The unique nonlinear optical effects in waveguides have recently become an 
area of intensive research (Ref. 24). Because of the optical confinement in a 
waveguide, guided optical waves can maintain high-power densities over a long 
distance for efficient nonlinear optical interactions at relatively low total 
power levels. Many important optical systems, such as optical communications 
and signal processing systems, involve guided-wave optics. Many desirable 
optical operations, such as optical PC, optical bistability, and optical 
switching, require nonlinear optical interactions. It is therefore important 
to understand the nonlinear optical wave behavior, including the propagation 
and coupling of the optical fields, in the waveguides. 

There has been some recent work on the propagation and coupling of optical 
waves in nonlinear waveguides (Refs. 24-26) and nonlinear couplers (Refs. 27 
and 28). Usually the techniques used are coupled-mode theory (Refs. 29-31) 
and beam-propagation method (Refs. 32-34). Strictly speaking, the concept of 
normal modes is applicable only for linear waveguides. The extension of the 
coupled-mode theory to nonlinear waveguides is a good approximation only when 
the nonlinear polarization is relatively small compared to the input fields. 

On the other hand, the beam-propagation method is not limited by this 
restriction and can be a more accurate approach for nonlinear waveguides. 
However, the beam-propagation method is formulated in terms of a scalar field 
(Refs. 32-34). This restricts the method to cases in which the electromag¬ 
netic (EM) fields can be separated into linearly polarized components, such as 
the transverse electric (TE) and transverse magnetic (TM) modes in slab 
waveguides, so that scalar-wave equations can be used to describe the problem. 
In many practical cases, such as optical-fiber waveguides with a circular or 
elliptical cross section and waveguides with a strong confinement, the 
polarization of the guided fields can become very complicated. Then the beam- 
propagation method becomes limited while the coupled-mode approach may still 
be good so long as the nonlinearity is not very large. It appears that each 
of these two methods has its own advantages and limitations and should be used 
in different situations. 
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Past work applying the coupled-mode formalism to nonlinear waveguides has 
assumed a real, isotropic dielectric constant and negligible contribution to 
the nonlinear polarization from the longitudinal field components 
(Refs. 24-27). However, many practical optical waveguides are fabricated with 
anisotropic materials, such as LiNb0 3 . Because of the tensorial nature of 
nonlinear susceptibilities, the longitudinal component of the nonlinear polar¬ 
ization can be significant in a nonlinear waveguide. Recently, there is also 
much interest in nonlinear waveguides with an active gain, i.e, with a complex 
dielectric constant, such as III-V semiconductor waveguides under current 
injection (Refs. 14, 17, and 35). This report considers the coupled-mode 
approach for nonlinear waveguides formally for various conditions of interest, 
including anisotropic waveguides and waveguides with a gain or loss. 

3.2 WAVE EQUATIONS 

The EM fields in a nonlinear waveguide are governed by the following Maxwell's 
equations, 


VxE - - i | (1) 

C dt 

VxH - i + (2) 

c dt c dt 

V • (E + 4*P) - 0 (3) 

V • H = 0 (*) 


where the material is assumed nonmagnetic so that B - H . The total 
polarization is given by the sum of the linear and nonlinear polarizations 
(Ref. 25) 

P(r,t) =P L (r,t) +P NL (r,t) 
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where 


P L (r, t) = J x (1) (r - x', t - t') * E(r', tOdr'dt' 


( 6 ) 


P^Cr, t) = J x <2) (r - r 3 , t - t 1 ;r - r 2 , t - t 2 ) 

: E(r 1 ,t 1 )B(r 2 ,t a )dx 1 dr a dt 1 dt 2 

m 

* fx t ' !1 (r-r 1 ,t-t 1 ;r-r ! ,t-c ! r-r J ,t-c 3 ) 
i E(r 1# t 1 )E(r 2 , t 2 )E (r 3 , t 3 )dr 1 dr 2 dr 3 dt 1 dt 2 dt 3 + ... ( 7 ) 

At a given frequency w, the electric and magnetic fields are 

E(r,t) = E(r, w)exp(-iut) + E*(r,ii>) exp(iwt) (8) 

H(r,t) = H(x, o>)exp(-i«t) + H* (r, w) exp (iwt) (9) 


and the complex polarization, P(r,u), is similarly defined. In case the 
linear and nonlinear electric susceptibilities are temporally nonlocal but 
spatially local functions, they depend on the light frequency but not on the 
wave propagation constants. Then, one can write the nonlinear Maxwell's 
equations in terms of E (r,w) , H(r,a>) , and P^Cr^ ) in the frequency 
domain, 


V x E - — H 

(10) 

c 

V x H - - —€ * E - l4 -* Q P" 1 * 

(ID 

c c 

V • (€ • E + 4xP nl ) = 0 

(12) 

V • H = 0 

(13) 
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The linear and nonlinear polarizations in the frequency domain are given by 


P L (r,w) - x a) (r, 0 >) • E(r, u) 


P NL (r,w) = x <2) (t,u =<•>! + w 2 ) : E (r, u 1 ) E (r, o> 2 ) 

+ x (3) (r, q = + u 2 + o 3 ) : E (r, to x ) E (r , a> 2 ) E (r , oi 3 ) 


The linear dielectric constant is in general a 3 x 3 tensor defined by 


£ (r,w) = 1 + 4n% (1) (r,w) 


so that 

E + 4jiP l = £ * E 

in the frequency domain. Here the local spatial dependence of x and e 
accounts for the spatial variation of material properties in the waveguide 
structure. 

In what follows, take the z-coordinate as the longitudinal direction of the 
waveguide. The fields can be separated into the transverse components in the 
cross-sectional xy-plane and the longitudinal components parallel to the 
z-direction of wave propagation in the waveguide. 


B = E T + 


H * H T + H z 


and the operator v can be written as v - V T + d/dz. Notice that the third and 
fourth Maxwell's equations (Eqs. 12 and 13) can be obtained by taking the 
divergence of the second and first Maxwell's equations (Eqs. 11 and 10), 
respectively, and do not provide more information than that obtained from the 
first two Maxwell's equations alone. Therefore, in most situations do not 
invoke them in the formulation of the problem. 
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3.3 NORMAL MODES OF THE LINEAR WAVEGUIDE 


Without the nonlinearity, there exists a complete set of normal modes of the 
waveguide structure. When e(r,u ) does not have z-dependence, at frequency u> 
each normal mode is characterized by a set of mode fields, 

E a (r,t) = f a (x,y)exp(iP a z - iwt) + (x-,y) exp (-ip*z + io>t) (20) 

H a (r, t) = <p a (x,y) exp (iP a z - iwt) + <p* (x,y) exp (-ip’z + iut) ( 21) 

where /9 a is the propagation constant of the ath mode. These normal mode 
fields propagate in the z-direction and satisfy the linear Maxwell's equa¬ 
tions, which are obtained by setting P*”* - 0 in Equations 1-4, with 
appropriate boundary conditions determined by the waveguide structure. 

The complex mode fields of the ath mode at the frequency w are 

E a {r,w) « i|r a (x,y) exp(ip a z) = i|i aT (x,y) exp (iP a z) 

+ i|r az (x,y)exp(ip a z) (22) 

H a (r,u) = <p a (x,y) exp(iP a z) = <p aT (x,y) exp (ip a z) 

+ <P a z (x,y) exp(ip a z) (23) 


where \|r aT - + V’ayY and + a2 - ^ a3 ,z. These complex-mode fields satisfy 

Equations 10-13 with P NL - 0. With the complex normal mode fields in 
Equations 22 and 23 and from the first two linear Maxwell's equations 
decomposed into longitudinal and transverse components, one has 




-T—— V T X <p aT - — € z 
€ zz 


”aT 


(24) 


€ t ' + € T ' = - ITT ( V T X «Paz + iPa 2 x <PaT ) 


(25) 
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(26) 


<P.Z - 777 V T x ♦ 


lb) 


aT 


<PaT = 177 (V T X taz + iPa 2 x *a T ) 


(27) 


In Equations 24-27, e - e x + e y + t z - c T + « z was used with the following 
definition 

•= - *<«xxE, * - e„E s ) (28) 

Cy 'E =<»<€y,E x - EyyEy * ey,E,) (29) 

€, -E = 2(€„E X ♦ €„Ey » e„E,) (80) 

Since the coupled-mode theory is based upon the linear normal modes of the 

waveguide, the properties of these normal modes will be examined before using 
them to formulate the problem involving nonlinearity in the waveguide. 

3.3.1 Backward-Propagating Modes 

For each normal mode, a backward-propagation mode can be constructed by 
time-reversal or space-inversion in the z-direction or both, depending on the 
property of the waveguide material, i.e., the property of c. Since only 
nonmagnetic (reciprocal) materials are considered in this report, the 
dielectric tensor is symmetric, « k j - tj k (Ref. 36), but it can be anisotropic 
as well as complex. In the following, the subscript -a will be used to 
represent the mode propagating in the backward direction with respect to the 
mode with the subscript a which propagates in the positive z-direction. 

Certain useful properties of the modes can be obtained by considering the 
relationship between the forward- and backward-propagating modes. 

CASE 1 . c is real, c(w) - e*(w), but it can be a general 3x3 tensor with 
nonvanishing real elements. This includes the case when the waveguide 
material is real but anisotropic. 
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In this case, time-reversal can be used, but not z-inversion, to construct the 
backward-propagating modes from the forward-propagating modes. By replacing 
time t with -t for time-reversal in the linear Maxwell's equations for the 
normal modes, the following relations are obtained 




K 


(31) 


<P-a - -<Pl 

and 

P-a = ‘PI 


(32) 


(33) 


For the propagating modes, /? a is real and there is a corresponding 
backward-propagating mode with /L a - -/3 a . For the evanescent modes, /9- a = - 

if}' is purely imaginary and there is no separate backward evanescent modes. 


CASE 2 . c(z) - c(-z) or e is cylindrically symmetric with respect to the 
z-axis. The tensor elements of e(u>) can be complex, i.e., the waveguide has 
gain or loss. 

In this case, the longitudinal z-axis of the waveguide is along one of the 
principal axes of the material, or the tensor £ can be transformed to a form 
of cylindrical symmetry 


e xx e xy 0 


€y X Cyy 0 
0 0 6 


zz / 


(34) 


Then, z-inversion can be used, but not time-reversal, to construct the 
backward-propagating modes from the forward-propagating modes by replacing z 
with -z in the linear Maxwell's equations for the normal modes. The following 
relations are obtained 


^-aT * ^aT ' ^-a* * ~^az 


(35) 
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*P-aT *PaT 


*P -az “ *Paz 


(36) 


and 


P-a = 'Pa 


(37) 


where ^ a , <Pa> and P a “ P a + i£a are all complex in general. 

CASE 3 . e is both real and cylindrically symmetrical with respect to the 
waveguide longitudinal z-axis, i.e. f e(u) - e*(w) and e(z) - e(-z). This 
includes the simplest situation when £ is real and isotropic. 

In this case, both time-reversal and z-inversion apply and the results in 
CASE 1 and CASE 2 can be combined to obtain the following relations 


♦•T = *aT = -aT = Y-aT * 


(38) 


'I'az = “♦« = “^-az = 'k-az * purely imaginary 


(39) 


<P aT " <PaT = -<P-ar ■ -<P-aT * * eal 


(40) 


<P az = -< Paz * <P- a , = -<p-az = purely imaginary 


(41) 


and 


P_ a = ~P a - real for propagating modes (42) 

P_ a = P a = purely imaginary for evanescent modes (43) 

CASE 4 . e (u>) w e*(u) and «(z) * e(-z). 

The material has a linear loss or gain and does not have a cylindrically- 
symmetrical optical property with respect to the waveguide z-direction. In 
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this case, both time-reversal and z-inversion do not apply and no simple 
relations can be obtained for the mode fields and propagation constants 
between the forward- and backward-propagating modes. 


3.3.2 Orthogonality and Normalization 


The linear normal modes of a waveguide form a complete orthonormal set. 
However, the orthogonality relations between the normal modes depend on the 
property of the waveguide material. The orthogonality and normalization 
conditions can be derived for each case discussed earlier using the results in 
Equations 31-43 and the reciprocity of the fields in a nonmagnetic medium 
where e is symmetric. The Lorentz reciprocity theorem states that at any 
point for which the fields are source-free there is the relation (Ref. 37) 


V' (EjXH 2 -E 2 X Hj) = 0 


(44) 


for any two sets of EM fields at the same frequency existing in the same 
linear, reciprocal medium. Taking E x - ^> 1 exp(i/3-^z) , H x - <p 1 exp(i^^z) , and E 2 
- V’2 ex P(i^2 z ) ’ " <P2 ex P(i^2 z ^ f° r t * ie fields of any two normal modes and 

integrating Equation 44 over the waveguide structure, one gets (Ref. 38) 

(Pr + P 2 ) //dxdyS • (f 1T x <p 2T - \|r 2T x <p 1T ) = 0 (45) 

in terms of the transverse-mode field components. 

CASE 1 . c(w) - c*(w), real but anisotropic. 

Taking 1 - a and 2 - -b for Equation 45 and using Equations 31-33, one gets 

(p a - Pb) //dxdy2 • (t aT X <p^ T ♦ i*Tfa T X <p aT ) = 0 ( 46 > 
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If both modes are propagating modes, and /3 b are both real and 0 b = /? b . The 
following orthogonality relation is obtained 

//dxdy2 • (i}r aT x <p bT + i|^ T x <p aT ) = ±26^ (47) 

The factor 2 in the above equation is introduced to give the modes the 
following normalization, 

J J dxdyg * Re(iJr aT x <p aT ) = j Jdxdyi • Re(i|r a x <p a ) = ±1 ( 48 > 

In Equations 47 and 48, the positive sign is used if mode a propagates in the 
positive z-direction and the negative sign corresponds to mode a propagating 
backward. Notice that 6& - 0 if a - -b in Equation 47. The backward- 
propagating mode is orthogonal to its corresponding forward-propagating mode 
in this case with anisotropic media. 

For the evanescent modes, /9 becomes purely imaginary. From Equation 46, if 
one or both of the modes are evanescent, the integral in Equation 47 vanishes 
even when a - b. This is because both terms in the integrand of Equation 47 
become purely imaginary and are complex conjugate of each other when a - b for 
evanescent modes. Therefore, the normalization condition for the evanescent 
modes is 

f J dxdy2 • ijr aT X <p aT - J JdxdyS • i|r a x <p^ = i < 49 > 

The integral in Equation 48 vanishes for evanescent modes since the evanescent 
modes cannot propagate and do not carry power (Ref. 37). 

CASE 2 . e(z) - e(-z), cylindrically symmetric but complex. 

Taking 1 - a and 2 - -b for Equation 45 and using Equations 35-37, one gets 

(Pa - P b ) //• WaT x <PbT + x <PaT> * 0 (50) 
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Taking 1 - a and 2 - b , one obtains 


(Pa + Pb } //cbcdyS • (<(r aT X <p bT - t bx x <p aT ) = 0 (51 > 

Multiplying Equations 50 and 51 by £ a + y3 b and /3 a - P h , respectively, and 
summing the results, the orthogonality relation obtained is 

IJ dxdyS • i|r aT x <p bT - 0, if p a * ± P b (52) 

If one takes mode a to propagate in the positive z-direction and uses 
Equations 35-37, one can normalize the mode fields to get 

/ / ‘ taT x <Pa T * -//<*Xdy2 * tar x <P-aT 

= f J dxdy2 • 4r_ aT x <p aT 

= -f f dxdy2 • ijf. aT x <p. aT = 1 (53) 

Since the mode fields ^ and ip are complex in this case, normalization to yield 

Equation 53 will involve complex normalization constants. Notice that 
different from the orthogonality relations in CASE 1, the forward and 
corresponding backward-propagating modes are not orthogonal in this case and 
the orthogonality and normalization relations in Equations 52 and 53 are in 
terms of $ aT and f bT , instead of $ aT and 

CASE 3 . e(«) - €*(«) and e(z) - c(-z). 

In this case, both f aI and f aT are real and Equations 47 and 52 are equivalent. 
Also, Equation 48 is equivalent to Equation 53. Then the orthogonality and 
normalization relations are the same as those in Equations 52 and 53, but with 

/ f dxdy2 • * aT x 9 bT = J J*dxdy2 • *r aT x <p; T (54 > 


Equation 49 applies for evanescent modes in this case, too. 
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CASE 4 . c(w) * c*(w) and c(z) * e(-z). 


In this case, there are no simple relations between the field components 
between the forward- and backward-propagating modes. The orthogonality 
relations can only be derived directly from Equation 45. 

The orthogonality relations derived in this section are in terms of the cross 
products of the transverse-mode fields. Orthogonality relations can also be 
obtained in terms of the dot products of the mode fields in CASES 2 and 3 and 
are given in Appendix A. 

3.4 NONLINEAR COUPLED-MODE EQUATIONS 

The nonlinear fields are solutions of the nonlinear Maxwell’s equations 
(Eqs. 10-13). Since Equations 12 and 13 are not independent of Equations 10 
and 11, only Equations 10 and 11 have to be used, which can be decomposed into 
longitudinal and transverse components, 


E = - -r-S— V. x H T - — c, • 

l(t>€ 1 € z 

M ZZ ^22 


E t - i* pf 

21 e 2Z 


(55) 


£ T • E T + € T • E z ~ 


Vj, X, H z ♦ 2 X I - 4«Pt L 
1W V z dz ' 


(56) 


H = -lV.xE. 

Z 10) 


(57) 


H 


T 


c 

i<i> 



x E„ 


+ 2 x 


dE T 1 
dz , 


3.4.1 Normal Mode Expansion 


(58) 


Strictly speaking, the linear normal modes form a complete basis only within 
the space of solutions for the linear Maxwell's equations. The solutions of 
the nonlinear Maxwell's equations in Equations 10-13 do not span the same 
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solution-space as that for the solutions of the linear Maxwell's equations and 
cannot be precisely represented as linear combinations of the linear normal 
modes. However, if the nonlinearity is not large so that 


|47tP NL |< < |e * E 


(59) 


the nonlinear polarization Pjjl may be treated as a perturbation to the linear 
polarization. Then, expansion of the fields in the nonlinear waveguide as 
linear combinations of the linear normal modes of the waveguide structure can 
be considered as a good approximation to represent the fields in a nonlinear 
waveguide. If Equation 59 is not satisfied, the coupled-mode approach fails 
for the nonlinear waveguide and beam-propagation method can be used if the 
fields are sufficiently linearly polarized so that the wave equations can be 
reduced to scalar-field equations. 

When expanding the EM fields in terms of linear combinations of normal modes, 
the same set of expansion coefficients are normally used for all the electric 
and magnetic-field components to keep the mode identity, 

E(r,co) = lA 4 (z) \|7 a (x,y) exp (iP B z) (60) 

a 

H(r,w) = Ea,(z) <p a (x,y) exp(ip a z) (61) 

a 

where £, sums over all forward and backward discrete guided propagating modes 
and integrates over the continuum of radiation modes. The evanescent modes 
should also be included in these expansions for completeness. However, since 
the evanescent modes do not propagate and are important only near sources or 
discontinuities of a waveguide, these fields will not be considered in the 
coupled-mode equations. 

The nonlinear waveguide fields are described by Equations 55-58, while the 
normal mode fields satisfy Equations 24-27. The field expansion in Equa¬ 
tions 60 and 61, with the same expansion coefficients for all field components 
to retain the mode identity, forces P 1 ^" to vanish in order for Equations 24 
and 55 to be satisfied simultaneously. However, in general, P 2 * 0 in 
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nonlinear waveguides, particularly when the waveguide has a strong confinement 
in which E z is significant or when the waveguide material is anisotropic in 
which P 2 can be induced by ELj.. When P z * 0, the linear mode identity 
cannot be completely retained in the expansion of the nonlinear-guided fields. 
The mode identity has to be relaxed and the coefficients for the six electric 
and magnetic-field components cannot all be the same. This is consistent with 
the observation that the solutions for the nonlinear system do not necessarily 
span the same solution-space as that for the linear system. However, there is 
no unique choice to relax the mode identity. In principle, all the expansion 
coefficients for the six field components can be different once the mode 
identity is not retained. It appears that the choice depends on how well the 
nonlinear fields can be approximated and whether the expansion coefficients 
can be easily solved. The linear mode field components are then used solely 
as a mathematically convenient function set to approximate the 
nonlinear-guided fields. 

One simple choice is to relax E 2 when P z * 0. By assuming a different set 
of expansion coefficients, say B a (z), for E z while maintaining the same set 
of expansion coefficients A a (z), for the other five field components, 

Equations 55-58 can be satisfied while simultaneously fulfilling the 
requirements of Equations 24-27. Therefore, when P z * 0, one can choose the 
following expansion, 


E t = £A a (z)\|r aT (x,y)exp(ip a z) (62) 

a 

E z = Ea,(z) ijr az (x,y) exp (iP a z) - — Pf* (63) 

• c zz 

H t - LA a (z) <p aT (x,y) exp(ip a z) (64) 

a 

H z - £A a (z) <p aT (x,y) exp (ip a z) (65) 

a 
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From the definition P® 1 * in Equation 15, it is clear that E z cannot be solved 
directly from the amplitude coefficients, A a (z) , but that Equations 62-65 must 
be solved by iteration. When - 0, Equations 62-65 reduce to those in 
Equations 60 and 61 which can be solved directly. 

3.4.2 Coupled-Mode Equations 


The fields in the nonlinear waveguide satisfy Equations 55-58. For 
generality, assume * 0 and substitute the field expansion given by 
Equations 62-65 in Equations 55-58. Then, using Equations 24-27, the 
following are obtained 


£ -T~ exp(ip a z)2 X <p aT 


i4n(o 

c 




( 66 ) 


dA P 

£ —exp(iP a z) 2 x i|r aT = 4 tiV t x — 


(67) 


Using these two equations and the orthogonality relations obtained in 
Subsection 3.3, the coupled-mode equations will be derived for the nonlinear 
waveguide in each case discussed earlier. 


CASE 1 . e(w) - e*(w), real but anisotropic. 

Taking the sum of the dot product of -^jwith Equation 66 and that of f^ T with 
Equation 67, integrating the resultant equation over space, and using the 
orthogonality relation for CASE 1 in Equation 47, one gets 


_ 


dz 


= 2xexp(-iP a z) J j dxdy 


* Y7 

|<PaT " V T X 


-NL 
* z 


... iw a* . mi®* _ 

+ —r~ vaT —-— ¥st 

c c€ zz 


« • «T * 


( 68 ) 
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Using vector identities and Equations 24-27, it can be shown that 

»nl . / -MT \ 

10 ) 


<p!t • V T x 


10 ) 




pNL + 
az * z 


ce. 




,NL 


aT 


+ V. 




V € zz 


*<PaT 


( 69 ) 


) 


Substituting Equation 69 into Equation 68 for the first term of the integrand, 
the following coupled-mode equation for a nonlinear waveguide is obtained 


dAa 

dz 


i 2 no) 


exp (-i P a z) J J dxdy P 1 * 1, ♦ 4r a 


pNL 

+ 2irexp (-ip a z) J J dxdy V T • —— x <p 


aT 


(70) 


Since every realistic waveguide is transversely bounded by vacuum in which 
P 1 ^ - 0, the last in integral Equation 70 can be shown to vanish by convert¬ 
ing it into a closed line integral at transverse infinity. Therefore, one 
gets 

“ ± exp(-iP a z) J J dxdyP WL • i|f a (71) 

where the plus sign is chosen if mode a is a forward-propagating mode in the 
positive z-direction and the minus sign is chosen if mode a is a backward- 
propagating mode. 

According to the normalization with Equation 48, the power in each guided 
propagating mode a is 

P.(z) - | A a (z) exp (i p a z) | 2 « ~ |A a (z) | 2 (72) 

were f) a is real for any propagation mode a in this case. 
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CASE 2 . e(z) - e(-z), cylindrically symmetric but complex. 
In this case, e T • P? - 0 and Equations 66 and 67 become 

? x <P«t " - 

dA P NL 

£ —i exp(ip a z) 2 x ijr aT = 4nV T x -i- 
a az e 22 


( 73 ) 


(74) 


/ # 

where $ a - /?„ + i/? a is a complex propagation constant. In this case, the 

forward-and corresponding backward-propagating modes are not independent. 

Using the orthogonality relations in Equations 52 and 53, one finds 

-^exp(ip a z) - ^exp(-iP a z) - ±±22. / f dxdyPf* * * aT < 75 > 

-^exp(ip a z) + -~£exp(-ip a z) - - J j'dxdypf' • i|f az < 76 ) 

Combining these two equations, one gets 

■ i^exp(-iP.z) //dxdy(pf • - p?*„) <”> 

2 - - QQa) eX p(i(3 a z) JJdxdyP" 1 - • *. < 78 7 

Notice the difference between Equation 70, where the components of f a are 
used, and Equation 77, where $ a has to be used, because of the difference in 
the orthogonality relations, Equations 47 and 52, respectively, for the two 
cases. 
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The total power in the guided field is 


■J2 |A a (z)A b (z)exp(i(p a + p b )z) Jf dxdy2 • Re (ty aT x <p bT ) 


Since the orthogonality relation and the normalization condition are given by 
Equations 52 and 53, respectively, in terms of faT and fbT* m °de power is not a 
well-defined concept here. However, when the gain or loss in the waveguide is 

small such that 0* < f) a , and f ±aT are approximately real functions 

according to Equations 38 and 40. Then the power in each guided mode can be 
approximated by 


P±a = -^-|A ±a (z) | 2 expU2p''z) 

ZTt 


CASE 3 . c(w) - e*(w) and e(z) - c(-z). 

In this case, following the same procedure in CASE 2, Equations 75 and 76 and 
Equations 77 and 78 are still obtained, but with real propagation constants f) a 
- -/L a for propagating modes. The power in each guided mode is given by 
Equation 72. 

CASE 4 . e(w) e*(w) and e(z) »< e(-z). 

In this case, a simple coupled-mode equation cannot be obtained because there 
are no simple relations between the mode field components and no simple 
orthogonality relations. 


In Subsection 3.4, the coupled-mode equations have been written for the 
amplitude coefficients, A a , using the electric-field pattern of the corre¬ 
sponding mode. This is the general practice used in the literature (Refs. 24, 
25, 27, and 29), where real and isotropic nonlinear waveguides with negligible 
have been considered. Because the mode orthogonality relations depend on 
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waveguide material properties, our equations for more complicated materials 
have different forms. Equation 71 (for the case of an anisotropic, real 
dielectric tensor), contains $*, while Equation 77 (for cylindrically symme¬ 
tric, complex dielectric tensors), contains f aT and t a2 . The coupled-mode 
equations can be written in a general form using the electric-field patterns 
of the counterpropagating modes 

- ± exp (-iP a z) J JdxdyP NL • ( 81 ) 

The plus sign is used when mode a propagates in the positive z-direction, 
while the minus sign is used when mode a propagates in the negative 
z-direction. This equation applies to all the cases where the orthogonality 
relations can be derived. 

Most of the commonly used waveguide materials, including nonlinear dielectric 
crystals such as LiNb0 3 and all III-V semiconductors such as GaAs, have aniso¬ 
tropic nonlinear optical properties. In this case, depending on the waveguide 
geometry with respect to the crystal optical axes, the TE-field components, 

E t , can induce a large longitudinal nonlinear polarization, P^. In addition, 
to increase the power density for efficient nonlinear interaction, a strong 
optical confinement is usually needed in nonlinear waveguides. The strong 
confinement can result in a large longitudinal electric-field component, E z , 
which can contribute to P z in the third-order process even if the material 
is isotropic. Therefore, in general, the longitudinal nonlinear polarization, 
Vf, does not vanish and can become very important in waveguides of special 
geometry. 

Subsection 3.4 shows that when P^ 0, the linear mode identity cannot be 
retained in the expansion of the nonlinear waveguide fields. In this case, 
there is no unique way of doing the expansion. In Subsection 3.4, the choice 
was to relax E 2 only. Another simple possibility is to choose one set of 
expansion coefficients, say A a (z), for all three electric-field components and 
a different set, B a (z), for all three magnetic-field components. By doing so, 
the coupled-mode equations become substantially more complicated and can be 
reduced to a simple form only under certain special conditions. This is 
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because the orthogonality relations, as given in Subsection 3.3, are all in 
terms of TE and TM mode fields, f aT and * aT , together. Therefore, using 
different expansion coefficients for Ej. and Ep makes it mathematically com¬ 
plicated. It is more convenient to choose the same coefficients for Ef and 
Hj.. In general terms, when the mode identity cannot be retained, it probably 
does not matter how the fields are expanded so long as the fields and the 
nonlinear polarization can be approximated as accurately as possible and the 
resulting equations are mathematically as simple as possible. It appears that 
the choice used in Subsection 3.4 does this better than other possible 
choices. If the material is magnetic with a nonvanishing longitudinal non¬ 
linear magnetization, then H z has to be relaxed and be expanded similar to E z 
in Equation 63. 

NL 

When P 2 - 0, the coupled-mode equations can be solved as coupled 
simultaneous equations for A a (z) directly since P x can be expressed explic¬ 
itly in terms of products of A a (z) through Equations 15 and 60. When P^" + 0, 
P 1 * 1 " cannot be explicitly expressed in terms of products or combinations of 
A a (z) any more because of Equation 63. Then, the coupled-mode equations 
cannot be reduced to simple and direct simultaneous equations for A a (z). In 
this case, the nonlinear fields can only be solved iteratively using the 
coupled-mode equations together with Equation 15 and Equations 62 and 63. 

A a ( z ) and pNL 

can be solved successively using the procedure outlined in 
Appendix B until a self-consistent solution is obtained. 

The contribution of Pf to the coupled-mode equation becomes important when 
Pf and f az are both significant, as can be seen from Equations 71, 77 and 78. 
This can happen when the waveguide material is anisotropic and the fields are 
tightly guided. When this contribution is negligible because P^* vanishes or 
because is very small such as in weakly guiding system, one has 

dA dA 

-^exp(iP a z) + —^exp{-i0 a z) = 0 < 82 > 

for CASES 2 and 3. In this case, taking a slowly varying amplitude 
approximation for the forward-propagating wave is equivalent to neglecting the 
backward-propagating wave (Refs. 25 and 40). This does not apply to CASE 1 
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since the forward and backward modes are independent in CASE 1. When the 
contribution P 2 is not negligible, Equation 82 is not valid. Then the 
relation between the forward- and backward-propagating waves can become more 
complicated. 

The linear and nonlinear electric susceptibilities, and the dielectric 
constant, e, are assumed all spatially local functions so that P L (r,u>) and 
P NL (r,w) can be written in the forms given by Equations 14 and 15, respec¬ 
tively. However, in certain cases of interest, such as in a semiconductor 
waveguide under carrier injection and the carriers have a diffusion length 
longer than the wavelength of the guided wave, % and € are nonlocal spatially 
as well as temporally. Then P L (r,w) and P^fr.w) have to be expressed as 
spatial convolution integrals of x and E. In this case, the z-dependent 
expansion coefficients, A a (z) , are also involved in the convolution integrals 
when the field expansion is substituted into the nonlinear wave equations. It 
seems that there is no simple coupled-mode approach in this situation. 
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4.0 EFFECT OF CARRIER DIFFUSION ON PHASE CONJUGATE FOUR-WAVE MIXING IN 
SEMICONDUCTOR WAVEGUIDES 

4.1 INTRODUCTION 

Four-wave mixing in a nonlinear optical waveguide has the advantage of long 
nonlinear interaction length and high field intensities due to the optical 
confinement of the waveguide. Phase conjugation by 4-wm in optical fibers was 
proposed by Yariv, et al. (Ref. 9). They showed that if a multimode optical 
fiber is excited with two counterpropagating single-mode pump beams both 
coupled into the lowest-order transverse mode, a third probe beam launched in 
the waveguide can be replicated with high fidelity through the PC 4-wm process 
in the waveguide. The multimode fiber has to accommodate a number of 
propagating modes equal to the number of resolution elements contained in the 
probe field for sufficient image resolution (Ref. 9). Subsequently, Hellwarth 
(Ref. 11) carried out a more detailed analysis on the generation of a PC 
replica of a free wave at the entrance plane of a multimode waveguide by 4-wm 
in the waveguide. The pump power required in the waveguide is orders of 
magnitude less than that needed for phase conjugation with unguided beams in a 
bulk medium. Furthermore, unlike for unguided waves, beam alignment is not 
critical (Ref. 11). It is only necessary to sufficiently couple the beams 
into the guided modes. 

Hellwarth’s analysis also showed that the two counterpropagating pump beams 
can be multimode without spoiling the fidelity of the PC replica, as long as 
the nonlinear coupling between the pump beams and the probe beam does not vary 
from mode to mode (Ref. 11). It is also possible to replicate the image 
carried by a multimode pump beam in the presence of a counterpropagating pump 
beam at the same frequency in the waveguide. In this two-beam configuration, 
the first pump serves as both the pump and the image beam. Both the three- 
beam and two-beam arrangements for phase conjugation were simultaneously 
demonstrated by Jensen and Hellwarth with a multimode glass waveguide filled 
with liquid CS 2 (Ref. 10). The pump powers required were nearly an order of 
magnitude lower than that required in the bulk medium. However, the conjugate 
reflectivity was less than unity because of the small optical nonlinearity of 
CS 2 . 
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Recently, there have been a series of experimental (Refs. 14, 17, 35, 39 and 
40) and theoretical (Refs. 16 and 41) studies of nearly degenerate 4-wm in 
semiconductor lasers and amplifiers. A large signal gain of 30 to 40 d£, 
which corresponds to a conjugate reflectivity of 10 3 to 10*, has been observed 
(Refs. 14 and 17). This is caused by the optical gain and the resonant 
enhancement in addition to the large optical nonlinearity in the semiconductor 
laser medium (Ref. 35). Using a noncollinear phase conjugation geometry in a 
broad-area diode laser, Lucente, et al. (Refs. 42 and 43) showed that the 4-wm 
signal is emitted in the PC direction. 

It becomes obvious that if the semiconductor waveguide can be used for the 
generation of optical PC signals, a very large PC reflectivity could be 
realized, in addition to the various advantages of using a dielectric wave¬ 
guide demonstrated by Hellwarth (Refs. 10 and 11). Since it is relatively 
easy to control the optical gain in a semiconductor waveguide by current 
injection, the signal intensity can also be easily varied without changing the 
pump intensity. So far, the experimental and theoretical investigations on 
the nearly degenerate 4-wm in semiconductor waveguides have focused on the 
nonlinear properties of the semiconductor medium without much consideration on 
the waveguiding and mode properties of the interaction (Refs. 14, 16, 17, 35 
and 39-42). Most of these studies have considered only single-mode semicon¬ 
ductor waveguides. However, as has been pointed out by Yariv, et al. (Ref. 9) 
and by Hellwarth (Ref. 11) for the generation of a PC image a multimode 
waveguide with a large number of guided modes is necessary for sufficient 
resolution of the image. In addition, the nonlinear coupling between modes 
has to be uniform for high-fidelity phase conjugation (Ref. 11). 

In a semiconductor waveguide under current injection for an optical gain, 
carrier diffusion can affect the efficiency of nonlinear wave coupling in a 
4-wm process (Refs. 41 and 43). In fact, using a noncollinear arrangement 
Lucente, et al. (Ref. 43) were able to measure the ambipolar carrier diffusion 
constant directly by measuring the spatial and frequency dependence of the 
conjugate signal efficiency in a broad-area laser diode. Carrier diffusion 
can also degrade the mode identity in a multimode waveguide. The effect of 
carrier diffusion on PC 4-wm in semiconductor waveguides is studied for the 
possibility and limitations of using semiconductor waveguides under current 
injection for efficient generation of PC image signals. 
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4.2 CARRIER-INDUCED POLARIZATIONS 


The 4-wm process is considered with two counterpropagating pump beams in a 
multimode semiconductor waveguide under current injection. The generic scheme 
is shown in Figure 2 where the z-direction is taken to be along the waveguide 
axis. The waveguide is assumed to be index-guided with a rectangular cross 
section of dimensions w and d in the x- and y-directions, respectively, as is 
shown in Figure 3. The intrinsic refractive index of the guiding core is n 
and that of the surrounding clad is n 0 < n for index guiding. The forward- 
and backward-propagating pump fields are Ep and Ep/, respectively. E i is the 
input image (probe) field and E B / is the returned signal (conjugate) field. 
Here a prime is used in the field indices for the backward-propagating fields, 
i.e., those propagating in the -z direction. The forward-propagating fields 
are indexed without a prime. For generality, consider nearly degenerate, but 
nondegenerate, 4-wm with degenerate pump frequency, w p - ay - w, and detuned 
probe frequency, a/* - w - fl. The returned signal is then at frequency 
ay - a> + Q. The degenerate case when o> p - ay - a> i - ay will be discussed in 
Subsection 4.4. 

The complex field at a particular frequency o>j is given by 

E.j(r,t) = Ej (r,u j )exp(-iw j t) (83) 

where Ej(r,«j) contains the spatial information, and the real field can be 
obtained from E.j(r,t) + E*j(r,t) . In a multimode waveguide, the spatially 
dependent parts of the fields can be expanded in terms of the normal mode 
fields of the waveguide. In a nonlinear waveguide where the longitudinal com- 
ponent of the nonlinear polarization, P z , does not vanish, the longitudinal 
z-component of the electric field can, in general, have different expansion 
coefficients from those for the transverse field components. In addition, the 
backward-propagating modes, f.,, are not generally identical to the forward- 
propagating modes, f a . However, for the discussions in this report, these 
differences will be ignored for simplicity in the formulation. Therefore, the 
spatially-dependent parts of the fields can be expanded as follows 

Ep(r) - (z)$,(x,y)exp(iP a z) (84) 
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Figure 2. Schematic diagram of PC 4-wm with a multimode waveguide. 
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Figure 3. Waveguide cross section. 
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( 85 ) 


Ep/(r) * E A< » > (z)t*<x,y)exp(-ip a z) 

a 

B^r) = (z) ^ a (x,y)exp( (ip a -fip a ) z) (86) 

a 

E s /(r) ■ (z ) ♦*(x,y)exp(-i (p a + 6p a ) z) * (87) 

a 

where is the propagation constant of mode a. Here, is assumed the 
normal mode of the intrinsic waveguide structure without current injection. 
Therefore, £ a is approximately real since the imaginary part, P a , contributed 
by the intrinsic waveguide loss is generally very small in comparison to real 
part, P a . In this case, the optical gain contributed by the injected current 
is not included in f} m . The contributions of the carrier-induced gain and 
optical nonlinearities are described by the total carrier-induced polarization 
P with linear and nonlinear components. The mode field pattern t a (x,y) 
depends only on the transverse coordinates and are characterized by two wave 
numbers, £ a and *j a , which describe the transverse spatial frequencies of ^ a in 
the x- and y-directions, respectively. For guided modes, these constants 
satisfy the following conditions 

*2 ♦ »i2 ♦ P2 - k 2 < 88 > 

and 

k 0 < P a < k (89) 

where k 0 - tiqu/c and k - nw/c. 

To obtain the carrier-induced polarization, one has to calculate the linear 
and nonlinear susceptibilities of the semiconductor medium in the presence of 
the injected current and the multimode optical fields. A formal and rigorous 
approach is to calculate the carrier distribution with the band structure of 
the semiconductor and then to employ the density-matrix formalism for the 
interaction of the optical fields and the medium (Ref. 44). This approach 
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will enable a more accurate accountability for the microscopic properties of 
the semiconductor medium, such as the band structure and the crystal symmetry. 
However, it is mathematically much more involved. Since the interest is 
primarily in the effect of carrier diffusion in this report, the microscopic 
nature of the semiconductor medium will not be considered. Instead, a pheno¬ 
menological approach will be taken with the rate-equation approximation. In 
this approach (Refs. 41-43), the carrier diffusion is described by an ambi- 
polar diffusion constant, D, and the intensity gain, g(N), as a function of 
the carrier density, N, as approximated by (Ref. 41). 

g(N) = a(N - N 0 ) (90) 

where a is the gain coefficient and N 0 is the minimum carrier density needed 
for optical transparency in the semiconductor medium. The total carrier- 
induced susceptibility can then be approximated by (Ref. 41). 

X (N) = - (b + i)g(N) (91) 

4710) 

where b is the antiguidance factor (Refs. 45 and 46). Since the interest is 
also in the possibility of a large signal gain, only the situation of small 
pump-probe detuning is considered (Ref. 41), | O | T S 1, where r is the carrier 

recombination lifetime. The 4-wm efficiency is much higher in this case, when 
compared to the case of a large detuning such that |Q|t >> 1 (Refs. 41-47), 
and the rate-equation approximation is valid. 

A procedure similar to that outlined by Agrawal (Ref. 41) is followed for a 
single-mode waveguide. However, because the 4-wm is being considered in a 
multimode waveguide, the mode properties (including intermode coupling) must 
be taken into account in the formulation. Begin by considering the spatial 
and temporal variations of the carrier density, N(r,t) described by the rate 
equation (Refs. 41 and 43) 
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dN(r, t) _ J _ N(r, t) 
dt ed -c 

+ DV 2 N(r, t) 


cn 


g(N(r, t)) |E(r, t) | 2 


(92) 


where J is the injection current density, e is the electronic charge, and D is 
the ambipolar diffusion constant. The total field in the waveguide is 


E(r,t) = Epexp (-iwt) + Ep/e>p(-iu)t) + E^xp(-i (to - Q) t) 
+ E 8 / exp (-i («+Q) t) 


(93) 


where E p ,E p /,E i and E e / are given by Equations 2-5, respectively. Taking the 
strong-pump approximation that |Ep ♦ Ep/| >> |EJ and |E S /| , one can write 

|E (r, t) | 2 “ |Ep + Ep/| 2 ♦ {[ (Ep + Ep,) * E' 

+ (Ep + Ep/) * • Eg/] exp(-iQt) + c.c.) (94) 

where the first term is time-independent but can be strongly spatially 
dependent, as are the time-dependent terms in the equation, because of the 
multimode mixing. To solve Equation 92, one writes (Ref. 41) 

N(r,t) =N a (r) + [N m (r) exp (-iflt) + c.c.] (95) 

where N s contains the effect of gain saturation due to the pump fields and N m 
accounts for the contributions from the wave mixing with the probe image and 
conjugate signal fields. Substituting Equations 94 and 95 in Equation 92, one 
obtains 


L 2 V 2 N e = (1 ♦ |Ep + Ep/| 2 /I s ) Ng - (N (1) + N 0 |Ep + Ep/| 2 /l 8 ) (96) 


L 2 V 2 N ra = (1 - iQt + |Ep + Ep/| 2 /l.)N m 




(N 8 - N 0 ) [ (Ep + Ep/) * El + (E; + E*/) • Eg/] /l 8 


(97) 



where L - y®r is Che carrier diffusion length and I s = 2ichci>/ (cnta) is the 
saturation intensity scaled to the unit of |E| 2 . The carrier density in the 
absence of the optical fields is N (1) - Jr/ed, which gives the unsaturated and 
unperturbed linear gain 

g (1) = a (N (1) - N 0 ) ( 98 > 

The corresponding unsaturated linear susceptibility is 

~<i) --£IL (b + i)g (1> (99) 

411(0 

In the presence of optical fields and carrier diffusion, Equations 96 and 97 
are usually approximated by spatially averaging the fields to account for the 
carrier diffusion and then by dropping the diffusion terms (Ref. 48). To more 
accurately trace the effect of carrier diffusion, this approach has not been 
taken. Instead, Equations 96 and 97 were solved by expanding N B (r) and 
N,,, (r) in spatially-dependent Fourier series containing terms of longitudinal 
and transverse spatial frequencies corresponding to those generated by mixing 
between field modes. The effect of carrier diffusion is to include in the 
wave-mixing efficiency a Lorentzian lineshape function 

2!(5 a ±S b ,Tl a ± Tl b ,p ft ± P b ;Q) = (l + [(£ tt ± £ b ) 2 

+ * V 2 + ( Pa ± L 2 

- iQt}' 1 (100) 


which depends on the spatial and temporal beat frequencies of the mixing 
modes. Therefore, each term in the Fourier series will have a different 
Lorentzian multiplier depending on the beat frequencies. With this 
understanding, the Fourier-series solutions can be written for N g and N m in 
the following concise forms, 


Ng = N u> _ (N„ (1) -N 0 )S£(Q)|E p ^ y[Vl s 

i + s*(o)£ a (i ( a p) i ( a p ' ) )/i fl 


(101) 
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N m = ~ 


- N 0 )S£<Q) [(Ep * y) • Ej+ (E^+E*/) • E,/] /I, 
1 + Sf(Q)r a (li p) + i‘ p,) )/i s 


( 102 ) 


where 2(0) - 2(0 - 0) and 


(1) + N 0 E a (I (p) ♦ I <p ) )/I t 
i + E a (i <p) + i (p/) ) /i. 


(103) 


is the background carrier saturation which does not depend on the mixing 
spatial frequencies of the modes. I <p> = |A ( | >) | 2 and I (p * = |A 1 | 2 are the 
spatial-frequency-independent average intensities in individual pump modes. 

The spatial-frequency dependencies of 2(0) and 2(0) vary with the mixing modes 
and can only be determined after the fields are expanded using Equations 
84-87. 

Using Equations 101-103 and Equations 90 and 91, one can write 


Xs * x ( »’ ~ 


X < a ) S£( 0 ) lEp+Ep/j 2 /!; 

1 + <£(0)E a <I <p) ♦ I 


(104) 


X. * - 


jffitftg) [(Bp * gp/) • Ei* 4 (g ; 4 e ;,) • a,,] /i, 
l + S£(Q)E.(I <P) + l <p 0 )/l a ) 


(105) 


where x* is the susceptibility corresponding to the saturated gain and x m is 
the susceptibility responsible for the wave mixing. The saturated carrier- 
induced background susceptibility is 


„<U _ 

X B “ 


i ♦ E a (i <p) ♦ i (p> )/i s 


(106) 
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In the strong-pump approximation, one can then write the total carrier-induced 
polarizations at the pump, image, and signal frequencies, 


P(o>) = P P (u) + P B (u) =x 8 (Ep + Ep,) 


(107) 


P(G) - Q) = P P (w - Q) + P b ( 0)-Q) = Xs ♦ X»(Ep + Ep/) 


(108) 


P(w+Q) = p p (0) + Q) + P B (« + Q) » X 8 E./ + X»(Sp + Ep/) 


(109) 


where P F and P B are polarizations containing forward- and backward- 
propagating wave vectors, respectively, i.e., 


P F - sum of exp(ipz) , exp(i3Pz) , . . .terms 


( 110 ) 


P B - sum of exp(-iPz) ,exp(-i3Pz) , . . .terms 


( 111 ) 


The total carrier-induced polarizations given in Equations 107-109 contain the 
linear and nonlinear polarizations. In general, the total pump intensity can 
easily be comparable to or larger than the saturation intensity, I s , in a 
semiconductor waveguide, which is typically of the order of 3 to 5 mW//m 2 
(Ref. 49). Therefore, the wave-mixing process has to be described by 
Equations 107-109 together with the total x s and Xm given in Equations 104 and 
105, respectively. One can expand x« and x m . 

x, - x'i’ * xT(<*p • K) ' ( e p- • O ♦ <*p • *;.) * • e;,)) <U 2) 


X. - xT ((*p • *;) * • b;) - (e; •«.,) - (e;, • e.,)) 


(113) 
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where 


v <3> _ 

J. 8 " 


_ a(Q)xi 1 ) /i s _ 

1 + S£(0)I a (li P) + li p/> )/l 8 


(114) 


v <3) .. 
a m “ 


a(Q)x < B ) /i e 

l + $£ (Q) E a (I (p> + I (p,) )/I. 


(115) 


The symbol for example, indicates that the spatial-frequency 

dependencies of £(G) in X*^* for the first term in Equation 113, 

depends on the spatial beat frequencies of the mode components 
of the multimode fields Ep and E* . It does not mean that spatial averaging 
has been taken for Ep * E^ . Using Equations 112 and 113, one can expand 
P (o>) ,P (o-Q) , and P(u + Q) in Equations 107-109, respectively. Since 
the phase-matched or nearly phase-matched nonlinear are the only interested 
processes, keep only terms -exp <iflz) for P p and those ~exp ( ' i ^ z) for P B and drop 
terms with higher-harmonic propagation constants. Then one has 


*r(»> » X l »\ + Xs 3, ((Ep ‘ E* p )Ep + (E p , * E^ ♦ (Ep • E^Ep,) ( 116 > 


P p (q - Q) « x'1% + x‘. 3) ((Sp * E*p)E i + (Ep, • E^Ei) 

♦ X m’ *((®*P • + (K> • E i)Ep' + ( E P • K/fy 

+ <v • e ;x) 


(117) 


P F (0) + Q) » Xs 3) (Ep * EJ,)b,, + X» 3> ((Ep • Ei)Ep + (e;, • E s ,)Ep) < 118 > 


».<«> = X^Ep, ♦ xi 3> ((Ep * KX>' + (Ep/ • E*/)Ep, + (E; • Ep,)Ep) 


(119) 


P B (w - Q) = x'. 3) (E‘p * Ep,) E d + X» > ‘((E*p • K i )E p , + (Ep, • E^,) < 120 > 
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(121) 


»,<*»♦ o> * x'. 11 * .< * x‘.”((Sp • - <*p' • e;X') 

* x» 3> «*„ • »;>v - <®p' • + 

*<k;, • Eg/^Ep/J 


Here P p (o)) and P b (q) contribute to the forward and backward pumps, 
respectively, P p (b> - Q) and P B (o> + Q) and contribute to the probe image 
and returned signal fields, respectively. The other two polarizations 
P f (g) + Q) and P B (w - Q) generate, in addition, a forward-propagating field 
at u + B and a backward-propagating field at u - ft through the 4-wm process. 

In the nondegenerate case, these two additional fields can be easily separated 
from the image and signal. In the degenerate case, they can become additional 
sources of complications for the PC process. This is discussed in 
Subsection 4.4. 


With these carrier-induced polarizations, the mode-amplitude coefficients of 
the multimode fields in Equations 84-87 satisfy the following coupled-mode 
equations. 


dA 1 _ j2nfa) 

dz C 


exp(-iP a z) JJ dxdyP F (w) • t|r a 


( 122 ) 


dA 


<p'> 


i2xw 


dz 


exp(iJJ a z) f f dxdyP B (u) * i|i a 


(123) 


= -~ 2 ”iSL ~ Q . } - exp(-i (P a -«P a )z) f JdxdyP F (o> - Q) - ijr a 


(124) 


—~— = - x2lt ( ^ - — exp(i (P a + 8p a ) z) f f dxdyP B (w - D) • t|r a (125) 

In these equations, the total carrier-induced polarizations, P F = P ( F ’ + P ( F ’ 
and P B = P^b* + P ( b > . which include the carrier-contributed saturated optical 
gain in P f p ' and P^* were used. This is because the defined waveguide modes 
were in terms of the intrinsic waveguide structure without injection, as has 
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been discussed in the previous section. In the strong-pump approximation, if 
the pump intensities are kept constant, the carrier-induced background 
saturated susceptibility, x's* > does not vary with the mixing fields. Also 
included is x's* in the total linear dielectric constant of the waveguide used 
to define the modes. Then P F and P B in each of Equations 122-125 should be 
replaced by P ( |* and P^* , respectively. However, /9 a and then become 
complex since the total linear dielectric constant becomes complex with the 
linear gain included. Both approaches give the same result, but the first one 
is more convenient for the analysis in this report. 


The nonlinear mixing terms in Equations 116-121 stem from various physical 


s. 

There are 

four 

(Ep 

•■&).« 

•v> 

K 



(V 

•*;><(«£ 

■ V) 

<Ep 

• *;>■(=»' 

• E v> 



•*:-) 

(Ep 


•«.') 


moving grating -cos (2 0z + fit) 


(126) 


(127) 


oscillating index - cos (Qt - 60z) 


(129) 


Notice that in addition to the z- and t-dependencies, every term has x- and 
y-dependencies which set up transverse standing gratings. For the ^Ep ' Ep) 
and (V • *v) terms, the background saturation contributions, which do not 
depend on the spatial frequencies, have already been accounted for in Xs X) and 
are included in P ( p' and P*^ . Therefore, these terms in P^' and P^' 
include only those transverse gratings set up by mixing between different 
modes in the pump fields. They disappear with single-mode pumps. The 
standing and moving gratings in Equations 126 and 127, respectively, are 
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created by interference between counterpropagating fields while the index- 
saturating transverse gratings and the time-dependent index oscillations, in 
Equations 128 and 129, respectively, are generated by copropagating fields. 

In the absence of carrier diffusion, such as 4-wm in an optical fiber, the 
spatial-frequency dependencies of the mixing efficiency through the Lorentzian 
function in each term disappear completely. In this case, every term in these 
equations contributes and the symbols, ^ and for the spatial-frequency 
dependence can be dropped. Where carrier diffusion is present, the spatial- 
frequency dependencies are important and the four categories of contributions 
have very different efficiencies. In a typical semiconductor waveguide, such 
as one of GaAs, n 0 - 3.5, A - 2irc/w -0.8 pm, and L - 2 pm. From Equation 89, 
this results in 20L > 2k 0 L > 100. 

Therefore, the nonlinear terms contributed by the standing and moving gratings 
vanish in the presence of carrier diffusion because of their vanishingly small 
efficiency according to Equation 100. Notice, that unlike the conclusion 
derived from spatial averaging discussed by Bogatov, et al. (Ref. 48) and that 
the effect of the moving gratings depends on the detuning frequency fi, the 
analysis shows that the moving gratings vanish irrespective of the detuning 
frequency, as can be seen from Equation 100. One conclusion can be immedi¬ 
ately drawn from these discussions: In the presence of carrier diffusion, the 
copropagating fields in the 4-wm process cannot be orthogonally polarized, 
while the relative direction between polarizations of counterpropagating 
fields does not affect the wave-mixing efficiency at all. For maximum 
efficiency, the copropagating fields have to be polarized in the same 
direction. 

Equations 124 and 125 describe the coupling between the probe image field and 
the conjugate signal. The relevant carrier-induced polarizations are 
P F (b) - Q) and P B (to + Q) with the standing and moving grating terms dropped 
from Equations 117 and 121, respectively. For the signal field (E 6 / in 
Eq. 87) to be the PC of the image field (E i in Eq. 86) with high fidelity, it 
is necessary that A a - K a A a for every mode a with a mode-independent 
constant, K, - K. The generation and development of the amplitude for each 
mode in the signal field is governed by 
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(s') 

~ - l2 * <“ exp(i(p. + sp,)z) j f dxdyi|r, • {xi 11 E s 

* x"’((«p -K) ♦ (Ep. •*;-))*.' * X» 3> (Ep * Ep» 

* xT<e;, -B.,)*,,} 

« - {SG + SM + PC + SC} 

c 


(130) 


The signal is first generated by scattering off the backward-propagating pump 
by the oscillating index fluctuation created by interference between E p and 
Ei. This process is described by the PC term. The first (SG) term provides 
the growth of the signal amplitude through the saturated gain, which is 
partially responsible for the very large signal efficiency in the 4-wm 
experiments with the semiconductor active medium (Refs. 14 and 17). The 
second and third terms (SM) account for the scattering and modulation of the 
signal induced by the transverse standing gratings created by the mixing 
between the modes of the strong forward- and backward-propagating pump fields, 
respectively. Once the signal is generated by the PC term, it is amplified 
and modulated by the SG and SM terms. In addition, it can beat with the 
backward-propagating pump field to create an oscillating index fluctuation 
which scatters off the backward pump field E p / into the signal. This process 
is described by the last (SC) term. The SG and SM processes do not create 
complexities in cross coupling between different modes in the image and signal 
fields, while the PC and SC processes are the source of complexities. The 
effects for these two terms are similar. To understand the effect of carrier 
diffusion, it suffices to examine only the PC term. 


With the multimode fields given in Equations 84-87, one has 


PC = f f dxdyXm 3} (Ep • Ej) Ep/ • i|r a exp (i (p a + 8 p a ) z) 

= C aaaa A a p ) A^ p) A^ i) *exp (i26P a z) + £ C abba A^ p ' ) A^ p, Ai i> *exp (i2&P a z) 

b#a 


+ £ C aabb A^ p ' ) A^ p> A ( b > *exp(i(8p a + 6p b )z) 

b#a 


(131) 
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where only the nearly phase-matched terms are kept. The subscript indices of 
the nonlinear coupling coefficients, and C aabb , correspond to the modes in 

the four mixing fields in the s'p'p i sequence. These nonlinear coupling 
coefficients are functions of spatial and temporal beat frequencies which 
depend strongly on the mixing modes, 


'abba 


„<1> 

Xa 


ff dxdyi|r a -i|r t 


se(P b - P. + &Pa*°> 


1 +2(P b -p a + SP a ;Q)£ a (li P) +li p/) )/I 


(132) 


'aabb 


„<1) 

X s 


/ f dxdyt|r a • ijr.j 


££(6P b ;Q)ilr b • 


1 + S£(»P b ;Q)E,Ui 


(p) 


♦I.* 0 )/!, 


(133) 


where Equation 115 was used and the transverse-spatial-frequency dependencies 
of £ are determined by the transverse patterns of the modes in the angle 
brackets. The dependence of £ on 5/3 in Equations 132 and 133 can be neglected 
since 60 L - nflL/c - 10~ 5 for flr - 1. C aaaa can be obtained by setting a - b in 
either Equation 132 or 133. Without loss of generality, it is also assumed 
that the injection current is uniformly distributed so that x's* is a 
spatially-independent constant in the guiding region. To illustrate the full 
spatial-frequency dependencies of the nonlinear coupling coefficients, take 
the mode pattern of the strongly index-guided mode (Ref. 50) as an example. 
With the origins of the transverse x- and y-coordinates defined at the center 
of the waveguide, the transverse mode pattern in the guiding regions, -w/2 < x 
< w/2 and -d/2 £ y < d/2, can be written (Ref. 50). 

¥ a (x,y) = -~cos (5 a x + a a ) cos (Tj a y + Y a )e a (134) 

v/wa 

where e a is the unit polarization vector of mode a, and w and d can be 
considered effective cross-sectional dimensions of the guiding region 
(Ref. 50). For simplicity, neglect the exponentially decaying fields outside 
the guiding region in the integration for C t bb , and C aij!)b . Since the interest 
is in multimode waveguides with many guided modes, most modes will be well 
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guided and this simplification will not change any general conclusions. It is 
also assumed that the copropagating fields are polarized in the same direction 
for maximum 4-wm efficiency. 

4.3.1 Low Saturation, Z a (li p) + li p ^/I, « 1 


In this case, since |g| £ 1,1 + S£(Q) £. (I ! P) + I a <p/) )/l B * 1. For a * b, 
one obtains 


'MU 


„<1) 

A a 

4wdl, 


{$£ (25 a , 2rj a , 0; Q) + 2$£ (2£ a , 0 , 0 ; Q) 


+ 2S£ (o, 2ri a , 0; Q) +4} 


(135) 


v (l) 

A B 


-abb. ' - 4 ^^ ( *a + + ’lb'P. ' Pb> 

+ 3 (5. ♦ 5 b ,n. - n b< P. - p b ;Q) (1 + 6^ ^) 

+ Sf (C. - S b ,n. + Tlb'P. ~ Pb' Q > ( 1 + 

- se <C., - « b .n. - n b ,P. ~ P b ;Q) d + *< a ,< b > d + 


(136) 


'aabb 


v (l) 

A 8 

4wdl, 


{2S£(0,2r lB/ 0;n)6 na>nb 


+ 2S£(2 5 a , 0, 0; Q)8 { . ({b ♦ 4) 


(137) 


From Equations 130 and 131, the mode intensity in the conjugate signal I^ s } a 
|A<*'>| 2 a|C| 2 ot|iff| 2 . From Equation 100, if any one of the following 
conditions, 

(i) |5 a ± 5 b |L £ 1. (ii) In. ± n b |L * 1- or (iii) |p a ± P b |L s l (138) 


is not valid, |iff| 2 will be reduced to <0.25. If any two of these three 
conditions are violated, then |iffJ 2 < 0.1. Using Equations 88 and 89, it can 
be shown that if Conditions (i) and (ii) are both satisfied, |/? a - 0 b |L < 1. 
Since the t functions in Equations 135-137 depend only on /9 a - 0 b and not on 


43 








/3 a + Pb> only Conditions (i) and (ii) need to be considered in the evaluation 
of the contribution of each term C aaaa , and C aabb to the conjugate signal. 

For the strongly index-guided modes given in Equation 134, one has - q a ?r/w 
and »j a - m a ir/d, where q a and m a are positive integers. The numbers of guided 
modes in the x- and y-directions are Q - 2nw/A and M - 2nd/A, respectively, 
and the total number of guided modes is N - QM. With n - 3.5, A -0.8 jxm, and 
L = 2 pm Conditions (i) and (ii) in Equation 138 become 

(i) |q a ±q b | * 0.018Q and (ii) |m a ± m b | s 0.018M (139) 

These conditions are hard to satisfy for the terms with spatial sum-frequency 
dependencies in Equations 135-137 unless Q and M are extremely large and a and 
b are modes of very low order. Therefore, except for a few low-order modes, 
one will expect 


'MM 


- C 


Mbb 


„(1) 

A B 

wdl„ 


(140) 


^-abba 


4wdl s ^ ~ ^ b ' 1,4 ” ^ b ' P* " Pb ' Q ) i 1 + 

C 1 + ^n..n b ) 


(141) 


Notice that for C .y. K . in Equations 131 and 141, the mode index a refers to the 
modes in the image and signal fields and the index b refers to the modes that 
exist simultaneously in the two pump fields. For the spatial difference- 
frequency dependencies in Equation 141, the conditions in Equation 139 can be 
satisfied only when a and b are neighboring modes. However, even when a and b 
are directly adjacent modes such that |q a - q b | - 1 or jm a - n^l - 1, Equa¬ 
tion 139 requires a minimum number of guided modes, > 55 or > 55, 

respectively. 

From Equation 131, the desirable PC signal 5 - K^ 15 * comes from the 
contributions of the first two terms associated with C aaaa and C^j,,, respect¬ 
ively. The third term associated with C aabb actually creates cross-talk 
between different modes in the image and signal fields and is not desirable, 
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term 


although it usually exists in the case of multimode pumps. The C aaaa 
accounts for the generation of the mode amplitude A‘ E } in the conjugate signal 
when both pumps are also coupled into the mode a. The term associated with 
C lbh . describes the more general source of generation of A| s 5 by all the other 
modes in the pump fields. In the case without carrier diffusion, this term is 
the major source of contribution (Ref. 11). Then, either single-mode (Ref. 9) 
or multimode (Ref. 11) pumps will work well to generate a PC signal. With 
single-mode pumps, the complication due to the C aabb term also disappears. 

These nice features are not available in the presence of carrier diffusion 
because of the requirement of Conditions (i) and (ii) in Equation 139 on C^,. 
As an example, assume that both pumps are coupled into only one mode, say the 
fundamental mode, and the waveguide supports a large number of guided modes in 
both the x- and y-directions such that N - QM » 55 2 - 3,000. Then, out of 
the total N available modes, the signal will be generated through C,^ only 
within the -3 x 10** N modes closest to the pump mode. With multimode pumps, 
the 0^, term only creates signal in modes around the pump modes. From 
Equation 140, C aaaa - constant and K a a A^ p } A^ p> , except for a few very low- 
order modes. Therefore, although the C aaaa term always exists, the signal 
amplitude A<* 1 in each mode depends on the amplitudes of the pump fields in 
the particular mode. From the discussions in this section, the conclusion is 
that the diffusion of carriers in a semiconductor waveguide severely compli¬ 
cates the coupling between various modes and limits the possibility of the 
generation of a high-fidelity PC signal. In fact, the only possibility of 
generating a high-fidelity PC signal with sufficient image resolution in a 
multimode semiconductor waveguide with carrier diffusion may be by coupling 
both pump fields into all the guided modes uniformly. Then the C aaaa and C^i,,, 
terms will generate A‘* ) - KA^ n * for all the guided modes a with K depending 
little on the mode index. However, even in this situation, the C aabb term 
exists and contributes to some undesirable cross-talk coupling between modes. 

4.3.2 High Saturation, 2 a (I^ p) + I a p * ) )/I, > 1 


In semiconductor laser waveguides, it is common to have significant gain 
saturation. When 12(0)Z a (I^ p> + I^ p ) )/I,\ > 1, the £ functions in Equa¬ 
tions 135-137 for C aaaa , C^,, and C aabb , respectively, have to be replaced by 
J?/(l + ZE a (I^ p) + I^ p ^/I,}, as can be realized from Equations 132 and 133. 
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This makes C 4444 , C^, ba> and C 44bb look very complicated, but it does not change 
much of the conclusions obtained in the case of low saturation. However, 
since now Z 4 (li p) + lj p '>)/I 4 > 1, it is possible that |*S 4 (l‘ p> + I‘ p '>)/I s | > 1 
even when Conditions (i) and (ii) are not satisfied and |jff| is small. The 
effect of gain saturation is to relax Conditions (i) and (ii) in Equation 138 
to 


<iv) |{. t J„|L s {e.(i * if' ) )/i.} 1/2 
(v) | n . i n b |L s {E.(i ‘f 1 * i ?'’) /1,} ,/2 U42) 


This is equivalent to the following relaxed conditions in place of 
Equation 139, 


(iv) jq a ± qj s 0 . 0180 {r a (l‘ p> + I < p '>) / I e } 1/2 
(V) |m a ± m b | s 0.0l8M{E a (l a <p) + I a (p,) ) / I 9 } 1/2 (143) 


Therefore, for those terms in C 4444 , C^, ba , and C 44bb which have spatial 
frequencies satisfying Conditions (iv) and (v) , Jf/{1 + JeE a (I^ p) + l^ p/) )/I s ) = 
l»/^«(l« p) + li p ) )• They become independent of the spatial beat frequencies. 
If the waveguide becomes so severely saturated that £ a (I a p) + I^ p 5 ) /I s » 1 
and Conditions (iv) and (v) are satisfied for the sum and difference beat 
frequencies between any two guided modes, the C^g and C aabb become totally 
independent of the spatial and temporal frequencies and, for a * b. 


’un 



4wdE a (li p) + 



(144) 


,U) 


'abba 


= c 


•abb 


4wd r, (i i p> * i?'>) 


(26, , ♦ 26. . 


4) 


(145) 


These are exactly what will be obtained in the case of no carrier diffusion. 
However, the nonlinear coupling coefficients and, therefore, the 4-wm 
efficiency are severely reduced by the increase of saturation. Notice that in 
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addition to the {£ a (I* p> + I^ p M)' 1 dependence in Equations 144 and 145, x a (1) 
is also inversely proportional to Z a (I a <p) + lj p ’ } ) according to Equation 108. 
Therefore, although gain saturation relaxes the stringent limitation of 
coupling the power of a pump mode to other modes in the signal due to carrier 
diffusion, this comes at a price of severely reducing the overall nonlinear 
coupling efficiency. In fact, \£\ falls off quickly when either of Conditions 

(i) or (ii) is not satisfied. For the transverse gratings generated by beat 
frequencies which satisfy Conditions (iv) and (v) but not Conditions (i) or 

(ii) , their contributions to the nonlinear coupling efficiency becomes unim¬ 
portant although Equations 144 and 145 apply. The conclusion is that gain 
saturation can relax the limitation caused by carrier diffusion but in a very 
undesirable way. For Equations 144 and 145 to apply for a large enough mode 
so that a high-fidelity PC signal can be generated with sufficient resolution, 
the saturation needed would be so high that the nonlinear coupling coeffi¬ 
cients become vanishingly small and the signal disappears altogether. 

4.4 DEGENERATE FOUR-WAVE MIXING 

In the degenerate case, when w p - w p » - - w,. , 0 - 0 and 6/9-0. Then, the 

moving gratings identified in Equation 127 become standing gratings like those 
in Equation 126 and the time-dependent oscillating index terms in Equation 128 
become only time-independent transverse gratings like those in Equation 129. 

In the presence of carrier diffusion, those terms with a longitudinal grating 
of spatial frequency 2/9 certainly do not contribute also in the degenerate 
case. Therefore, still only those terms identified in Equations 128 and 129 
are kept for the carrier-induced polarizations. 

As discussed earlier and expressed in Equations 124 and 125, in the 
nondegenerate case only P p (w-Q) and P b (u+Q) contribute to the generation, 
growth, and modulation of A^ n and A^*’respectively. The polarizations 
P p (o>+Q) and P b ((i>-Q) are also generated in the 4-wm process, but they can 
be easily distinguished from the contributions of P p (o)-Q)and P b ((i)+Q) and 
do not affect A^ u and A^ s \ directly. The major difference between the 
degenerate and nondegenerate cases is that in the degenerate case this 
distinction disappears and the coupling between A^ 15 and A^ s > is complicated 
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by these two additional polarizations, P p (b)+Q) and P B (w-Q) , as Cl-'O . In 
the strong-pump approximation and in the presence of carrier diffusion, one 
can then write 


P ( f P) («) = Xs% + xi 3> «Ep * + <Sp/ * E;/)Ep) (146) 

Pf 1 ’ («) = ♦ X ( « 3) ((Ep ‘ Bp)®i + (Bp' * Bp/^Ei) 

* X ( » 3) K ‘*0^ + (K -E^Ep) 

+ Xm’KBp ‘ ♦ (Eg, • E 6 ,)Ep) (147) 


Pa 15 '’ («) - Xs’V + xi 3, ((Bp • E;,)Ep, + (Ep/ * E;,)Ep/) (l48) 

Pb ' 1 (<*>) = Xa’B./ + X* 3> ((Bp * e;)B., + (Ep/ • E*/^E S /) 

+ Xm 3, ((Ep * El)Ep/ + (ig, • E a /)Ep/) 

♦ X* )m ({K * E,)Ep, ♦ (Ep/ • E*/^Ep/) < 149 > 


Equations 146 and 148 are identical to Equations 116 and 119, respectively, 
except that the terms corresponding to the standing longitudinal gratings are 
dropped due to carrier diffusion. Equations 122 and 123 are still valid with 
P F (u) and P B (u) replaced by P** 1 and P^* (w) , respectively. Therefore, in 
the strong-pump approximation, the pumps are not changed by frequency degen¬ 
eracy. However, the coupled-mode equations for A^ 15 and A* s) given in 
Equations 124 and 125 have to be modified by replacing P p ((«)-Q) and 
P B (<t>+Q) with P^* and P ( £ *(<*), respectively, in addition to setting fl - 0 
and 60 t - 0. Therefore, the generation, development, modulation, and coupling 
of A ^*’ ) is now described by 


(0 0 

—~ — * - i2n . P {SG + SM + PC ♦ SC + AM) (150) 

dz c 

where the SG, SH, PC, and SC terms are the same as those defined in 
Equation 130. 
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They also provide the same contributions and limitations to the PC 4-wm 
process as described in Subsection 4.3 for the nondegenerate case. The 
existence of the anomalous mixing (AM) terms in Equation 150 is due to the 
fact that P b (q-Q) given in Equation 120 also contributes to in the 

degenerate case, 

AM = ff dxdyXm ) * ((K * Ei) Ep/ + • E*/) E,,) • i|r a exp UP a z) (151) 

The effect of carrier diffusion on the multimode-field mixing of the AM terms 
in Equation 151 is similar to that discussed in Subsection 4.3 for the PC 
term. 

From Equations 150 and 151, it is obvious that the effect of frequency 
degeneracy is to couple A^ s#) to the modes of E i and E^/ directly. As a 
consequence, one can expect that A^*’ 5 becomes a sum of many different contri¬ 
butions. In addition to the desirable PC contribution, KjjA^*, discussed in 
Subsection 4.3, there is a non-PC contribution, K^A 1 * 1 + E^K^A ^ , to the 
backward-propagating signal amplitude A^*’There is also a backward- 
propagating amplitude proportional to A^*’due to the second AM term in 
Equation 151. Obviously, these extra contributions in the degenerate case 
tend to destroy the PC signal. With a single-mode waveguide, it can be seen 
that the returned signal in a degenerate 4-wm process can never be the PC of 
the probe field because the AM terms are as effective as the PC term in 
generating the backward-propagating signal in this case. However, as is 
discussed earlier, in order to have sufficient image resolution in a PC image 
replication process, a multimode waveguide with a large number of guided modes 
is necessary. It can be shown through a more detailed study of the AM term 
that the non-PC backward-propagating signals proportional to E A and E* a , are 
only generated in the modes where both pumps exist. Therefore, if the 
multimode waveguide is pumped with single-mode pump fields, all the modes in 
the signal will not be affected by the AM terms except the one mode where the 
pumps exist. However, this is only possible in a waveguide without carrier 
diffusion since carrier diffusion makes the coupling of the power from a pump 
mode to a different mode in the signal very inefficient, as has been discussed 
in Subsection 4.3. 
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On the other hand, it seems that even when the waveguide is pumped with 
multimode pump fields, the effect of the AM contributions is not to destroy 
the PC replica of the image but only to reduce the efficiency by coupling pump 
powers to a non-PC background noise signal. This is because only the PC 
components of the backward-propagating signal field can add together at a 
certain distance away from the waveguide entrance to reproduce the image, 
while the non-PC components, not having the correct total z-dependent phase 
relationship among the modes, tend to diverge and average out as a background 
noise at a distance. This probably explains why Jensen and Hellwarth 
(Ref. 10) could observe PC image replication through a degenerate 4-wm process 
in a waveguide pumped with multimode pump fields. 

Finally, there is a question as to how small 0 should be when it can be 
considered degenerate. To our understanding, so long as there are externally 
applied strong pumps as the arrangement shown in Figure 2, (1 is limited only 
by the ability to resolve the frequencies. In other words, the frequency 
degeneracy is instrument-defined. However, with a semiconductor laser wave¬ 
guide, it is possible to drive the laser above threshold and use the internal 
lasing fields as pumps. In this case, Q has to be kept larger than the 
injection-locking range (Ref. 17). Otherwise, the laser waveguide can become 
injection-locked and the PC process no longer exists. 

4.5 CONCLUSIONS 

Four-wave mixing in a semiconductor waveguide under current injection has a 
very large signal gain due to the large carrier-induced optical gain and 
optical nonlinearity in the semiconductor medium. The interest is in the 
possibility of utilizing this large conjugate signal gain for optical PC image 
replication. For the generation of a PC replica image with a multimode wave¬ 
guide, the participation of a large number of guided propagating modes is 
needed for sufficient image resolution. In addition, the PC efficiency has to 
be uniform across the modes in the image field for high-fidelity PC. Carrier 
diffusion in a semiconductor waveguide under current sets many limitations on 
these requirements. 

The effect of carrier diffusion on PC 4-wm in multimode semiconductor 
waveguides has been studied. Several conclusions are obtained from the 
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analysis. First of all, the gratings set up by counterpropagating fields in a 
4-wm arrangement are always washed out by carrier diffusion. The whole 
process relies on the transverse gratings and the time-dependent index oscil¬ 
lation set up by beating between copropagating fields to provide any conjugate 
signal reflectivity. Therefore, the copropagating fields cannot be orthogon¬ 
ally polarized. For maximum efficiency, they have to be polarized in 
parallel. On the other hand, the relative direction between polarizations of 
the counterpropagating fields does not affect the efficiency. The major 
effect of carrier diffusion is to severely limit the efficiency of coupling 
the power in a particular pump mode to other modes in the signal and image 
fields. In a typical semiconductor waveguide, the power in each pump mode can 
be coupled to only a very small fraction of all guided modes. Even in a wave¬ 
guide with a very large number of guided propagating modes, only modes very 
closely neighboring a pump mode can be generated in the conjugate field. This 
makes it totally impossible to pump with single-mode pumps. The only possi¬ 
bility is to couple both pump fields into all the guided modes with uniform 
amplitudes across the pump modes, which is certainly very difficult if not 
impossible. Gain saturation by strong pumps, which is common in semiconductor 
waveguides, relaxes this limitation. The power in a pump mode can be coupled 
to more modes in the signal field as the saturation increases. In the limit 
of very high saturation, the effect of carrier diffusion becomes unimportant 
and the nonlinear mode-coupling coefficients reduce to those of no carrier 
diffusion. However, this does not solve the problem since gain saturation 
severely reduces the 4-wm efficiency. To the extent that saturation suffi¬ 
ciently removes the effect of carrier diffusion, the conjugate signal also 
disappears altogether. 

The experimental observation of Lucente, et al. (Ref. 43) is consistent with 
our conclusions. In their experiment, PC 4-wm was performed with a noncol- 
linear geometry in a broad-area semiconductor laser waveguide. The PC signal 
was observed to have spatial- and frequency-dependencies of a Lorentzian 
lineshape function. In this experimental setup, the pump and probe fields 
were obviously coupled into only one spatial mode but were separated by the 
noncollinear geometry. The grating responsible for the generation of the 
conjugate signal was created by the noncollinear wave vectors of the forward- 
propagating pump and the probe fields. This is equivalent to the last term in 
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Among all the terms discussed in Subsection 4.3, only this term existed 
in their experiment because the noncollinear pump and probe fields were 
coupled into only one mode. 
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5.0 CONJUGATE GENERATION REQUIREMENTS IN MULTIMODE WAVEGUIDES AND THE 
UTILITY OF WAVEGUIDE ARRAYS 

5.1 INTRODUCTION 

The formalism developed in Sections 3.0 and 4.0 shows that some stringent 
requirements must be satisfied in order to generate a PC beam in a multimode 
waveguide. In the first part of the section, these requirements are discussed 
and placed in the context of previously published work. The additional 
constraints imposed when the nonlinear interaction is subject to diffusion 
effects are also considered. Finally, the idea of using 4-wm in an array of 
single mode waveguides as a means of overcoming the limitations imposed on 
multimode waveguides is introduced. The use of the array also allows optical 
control of the propagation direction of the conjugate wavefront, a feature not 
possible in multimode waveguides and only achievable in some bulk media using 
precise wavelength control of the interacting beams. 

5.2 MULTIMODE WAVEGUIDE PHASE-CONJUGATION REQUIREMENTS 

The first requirement for high fidelity phase conjugation is that the pump 
beams must not be depleted by the nonlinear interaction. As with bulk media, 
the "undepleted pump approximation" is necessary to keep the phase of the 
generated beam in each mode (or wave vector in the bulk case) from phase and 
amplitude variations due to the z-dependent amplitude coefficients of the pump 
beams. This is quite reasonable given the similarity between first-order 
differential equations for the coupled-wave and coupled-mode formalisms. 

Except for the overall proportionality factors and a replacement of the 
k-vector delta function in the bulk case with the transverse overlap integral 
of the mode patterns in the waveguide (Ref. 25), the two formalisms similarly 
relate the generated field amplitude to the input field amplitudes. 

A second constraint that generally arises relates to the fact the P 1 ^ must be 
treated differently from P T in the coupled-mode formalism. If the longitud¬ 
inal nonlinear polarization plays a significant role, then the coupled-mode 
equations must be solved iteratively, as shown in Section 3.0. This iterative 
process will not retain the phase relationship between the conjugate of the 
input signal and the generated beam and, therefore, lead to a degrading of the 
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phase conjugation fidelity. Though there are specific cases where P 1 ^ will 
not play a role due to the specific geometry or material properties, in 
general, phase conjugation in multimode waveguides is limited to those 
situations where it can be ignored. This requires that all modes involved in 
the nonlinear interaction be guided and that the condition of weak guiding 
hold (Ref. 29). Weak guiding requires that the modes vary slowly, with 
respect to the optical wavelength, in the transverse directions. Because the 
longitudinal component of the electric field is related to the transverse mode 
derivatives in Maxwell's equation, the z-components of the electric field and 
nonlinear polarization can generally be ignored when all modes involved in the 
interaction are weakly guided. 

Assuming that P z can be neglected, the mode expansion for the electric field 
leads to overlap integrals, 


C abcd = //dxdyi|r_ a * x (3) : i|f b ilr c ilr d (152 

and phase factors, expl-i(^.-^ b -/9 e -^ d )z), in the coupled-mode equations. In 

general, Equation 152 implies coupling among all combinations of modes in the 

mixing process. However, Hellwarth has pointed out that only pairwise mode 

combinations, e.g., (a,b,c,d) - (a,b,-a,-b) or (a, a,b,-b), etc., make 

significant contributions in the case where the waveguide length, i, is 

sufficiently large and the waveguide modes are nondegenerate (Ref. 11). In 

this case, the overlap integrals essentially become intensity overlap 

integrals plus a term which describes polarization overlap of the modes and 

C«bcd “ C«b. With a proper choice of input field polarization, all excited 

modes can be made collinearly polarized in the weakly guiding waveguide 

NL 

required for negligible P z (Ref. 29). 

The coupled-mode equation for the PC can be written in the form, 


^ - a — = - -i£!Ll5U_9l {SM+ PC} 
dz c 


(153) 
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where 


SM = 2C m (|A i p) | 2 + |A i p/) | 2 )a i“ 7) + 
+ (a b P> A a p> * ♦ A l p/) A { .i ] ’)A 


rcJlA^p ♦ | A i p/) | 2 )a i* 7 ’ 

b*a 

-b /} exp (i (&P a -&P b ) z)} 


( 154 ) 


PC = C aa Ai P) Al p/, Ai 1, *exp(i28p a Z) 

+ Ec ab {Ab P> A ( ^ ,) A a (i> *exp(i26p a z) 

b*a 1 

+ A^’A^A^^expCitfipa + 8P b )z)} (155) 

6/9 a and 6/9 b in Equations 154 and 155 arise from the frequency detuning, ft. 
Because intensity overlap integrals are not expected to show a strong mode-to- 
mode variation (Ref. 11), the first two terms in both SM and PC have the same 
form as the bulk case. The last term in SM and PC arise from the spatial 
overlap of pump and probe beams. Therefore, as in bulk media, phase 
conjugation by 4-wm in waveguides requires undepleted pumps and is only 
possible when can be ignored, i.e., the weak-guiding or slow transverse 

variation limit. If the image and pump beams occupy the same modes, then the 
image beam cannot be depleted in the 4-wm interaction. This is the case 
analyzed by Hellwarth (Ref. 11). Yariv, et al., analyzed the case when the 
pump and probe beams occupy different modes and showed that image depletion 
and conjugate amplification may occur (Ref. 9). Phase conjugation spatial and 
temporal bandwidth in a multimode waveguide are limited by the requirements of 
sufficient guiding so that only pairwise spatial overlap integrals contribute, 
and no dephasing of the spread of pump and image frequency components over the 
interaction length. These two requirements can be summarized: 

^ f >1 (156) 

S 


and 


2Q 

c 


{ < 1 


(157) 
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where S is the effective cross section of the waveguide. It can be approxi¬ 
mated (Ref. 29), 


S - 2n An A 


( 158 ) 


where A is the core cross-sectional area and An is the index mismatch between 
core and cladding. Putting Equations 156 and 157 together limits the number 
of modes to on the order of 10 3 at a gigahertz bandwidth. Jensen and 
Hellwarth observed phase conjugation in multimode waveguide containing on the 
order of 10 5 modes but the conjugate images they reported could have been 
formed from only the lower order modes of their waveguide (Ref. 10). 

In addition to the above considerations, diffusion must be considered in those 
situations where the nonlinearity results from mobile particles, such as the 
free carriers in a semiconductor under current injection. In a rectangular 
waveguide, each transverse mode field pattern, ^(x,y), is characterized by two 
transverse wave numbers, and rj a , associated with the transverse spatial 
frequencies of in the x- and y-directions, respectively. These wave 
numbers satisfy the following condition (Ref. 50), 


il * nJ * P2 - <: 

Due to the carrier diffusion, the nonlinear susceptibilities depend on the 
spatial and temporal beat frequencies through a Lorentzian function, 


± Sb.Tl, * n b ' Pa ± Pb'* 

= {1 + [<5. ± 5 b ) 2 + (11a ± 11b) 2 + (Pa * P b ) 2 ] L 2 - iQt}' 1 

where r is the carrier lifetime and L is the carrier diffusion length. 


As was discussed in Section 4.0, the contribution of a term in the nonlinear 
susceptibility becomes negligible unless the transverse beat frequencies in 
the t function of this particular term satisfy (i)j£ a ± £ b l L - 1 and 
(ii)|n« ± *} b |L - 1- These conditions set stringent limitations on the 
efficiency of coupling the power from a mode in the pump fields to other modes 
in the signal. As an example, in a typical semiconductor waveguide with 
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n - 3.5, A - 2irc/w - 0.8 pm, and L - 2 pm, each pump mode can only couple to 
the nearest neighboring modes within <0.024/y/An of the total guided modes in 
both x- and y-directions. This means that in a truly two-dimensional 
waveguide with a sufficiently large number of modes for image resolution, each 
pump mode will couple to only <6 x 10'VAn of the total two-dimensional guided 
modes. As a result, to uniformly couple all the guided modes, the index step, 
An, cannot be >6 x 10' 4 . This results in a very weakly guiding waveguide with 
a small acceptance angle of -7 deg at the input from the free space and a very 
small number of guided modes in the waveguide. When the index step is 
increased to improve the pixel resolution and the acceptance angle, the 
multimode coupling becomes highly nonuniform in a waveguide with carrier 
diffusion. 

5.3 FOUR-WAVE MIXING IN A WAVEGUIDE ARRAY 

To avoid the problems involved in a multimode waveguide while still keeping 
the advantages of a waveguide, 4-wm in an array of single-mode waveguides is 
proposed. The geometry of such an arrangement is shown in Figure 4. 



Figure 4. Geometry for 4-wm in a waveguide array with incident pump beams 
and E 2 , signal beam E s , and generated beam E c . 


A periodic array of identical waveguides is located with faces in the x-y 
plane at z - 0 and z - -L and with confined propagation along the z axis. 

Each waveguide supports only one mode with a propagation constant £, which may 
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be complex, and an electric-field pattern with the transverse component VKx.y) 
polarized along the y axis. Three beams are incident upon the array; two on 
the z - 0 face and one on the opposite face, with the constraints that their 
propagation directions fall within the acceptance angle of the waveguide and 
that their amplitudes are slowly varying across the array face with respect to 
the interwaveguide spacing. Two of the beams, one on each face, are pump 
beams at the same frequency, u 0 , while the third is shifted by an amount 0. 
Within the waveguide mode all beams are spatially degenerate, but 4-wm 
involving the probe beam is frequency shifted from the nonlinear mixing of the 
two pump beams. 

Coupling between the guided mode of each waveguide and the radiation fields is 
completely determined by the continuity of the tangential electric-field com¬ 
ponents at the interface and the constraints on the propagation constants in 
the guided region and free space. If all three incident fields have their 
tangential polarization along the y axis, to match the guided-mode polariza¬ 
tion, then the coupling of the three fields into the waveguides will be 
determined by their angles of incidence. The angle of incidence will also 
determine the phase relation between the fields incident at each waveguide. 

The three input waves are decomposed into a linear superposition of plane 
waves, 


■u (x,y, 2 ) = dk jk dk jy exp [ikj • x] £] (k jx ,k jy ) (161 > 

where the subscript j is s, 1 or 2 for the signal and first and second pump 
waves, respectively. In free space, 






(162) 


and 


- o 


(163) 
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allowing the z components of ^ and kj to be determined in terms of the x and 
y components. If the surface reflection can be ignored because of, for 
example, an antireflection coating, then the electric field from the signal 
beam coupled into the waveguide array at the interface will be (Ref. 51), 


E s (x,y, 0) = £ E, mn (x < y / 0) 

m, n 

- £ A smn (0)ijf (x-mx 0 ,y-ny 0 ) 

m,n 


(164) 


where x 0 and y 0 are the center-to-center element spacing in the array and m 
and n denote the individual waveguides. The amplitude coefficients are deter¬ 
mined by the mode-plane wave overlap, 


<' 0 > = //^sx^sy^sy <k«c, k sy> 

x ffdxdyty (x-mx 0 ,y-ny 0 ) 
x exp [i (k 8X x +k, y y) ] 

= / J^sx^syKy < Kx' Ky > 

x exp [i (mk BX x 0 +nk ay y 0 )] 

y(k xy ,k By ) = J Jdxdyijr (x,y) exp [i (k 8X x+k sy y) ] 

where is the tangential component of the input plane wave. Similar 
expressions hold for the two pump waves. Note that the phase of each plane 
wave at all waveguides is explicitly retained and that waves ranging over the 
waveguide acceptance angle are all coupled into the mode. Propagation of the 
coupled light within the waveguide is determined by the propagation constant, 
P, and the nonlinear interaction. All light not coupled into a waveguide mode 
is assumed to be absorbed. 

The usual approximations can be made for a nearly degenerate 4-wm interaction, 
slowly varying amplitudes and undepleted pump beams (Ref. 1). In the wave¬ 
guide, the nonlinear susceptibility * <3> is weak enough that the linear mode 


(165) 

(166) 

(167) 
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pattern ^(x,y) describes the TE-field patterns. The frequency range, Q, is 
small enough so that both x <3) and yf> are constant over the interval u 0 ± fi. 

The A-wm interaction for a single-mode waveguide has been solved to show that 
the generated field E c is proportional to the conjugate of the input field 
(Ref. Al). If there is negligible coupling between waveguides in the arrays, 
then 

E cmn( x 'y'°> = (0)i|f(x - mx 0 , y - ny 0 ) (168) 

Acrnn (0) = F (A^, A 2mn , 0) 0 , Q , P , L, D A^* (0) 

= FmnA 8mn *(0) (169) 

Field confinement effects are contained in the parameter T, 

r(Q) fdxdy x <3) (x,y, w 0 ,Q) |4»(x,y) | 4 (170) 

To determine the generated field propagating in free space, the previous 
procedure was followed, expanding the output mode patterns as a superposition 
of plane waves (Ref. 51). 


E c (x,y,z) -JJdk^dk^ 

x exp(ik c • r)£ ^(k cx ,k^) 

m,n 

^cnm^cx*^) = // E cmn (x f y,0) 

x exp [-i (k^pc + k^y) ] dxdy 


(171) 


(172) 


60 




Insertion of Equations 168 and 169 into Equations 171 and 172 yields, 

E c (x,y,z) = — ff dk^dkpy exp(ik c • r) 

4lt 2 JJ 

*ff dksx^sy Y* <k ax ,k 8y ) Y(-k cx , -k^) 
x K *<k flX ,k 8y ) £ F n« exp[-i(mk 8X x 0 

m,n 

+ nk 8y y 0 + mk^o + rak^)] (173) 

If is a constant, then the summation over waveguides reduces to a set of 
angularly separated 6 functions, corresponding to the diffraction orders of 
the array, when the array sufficiently samples the plane-wave components of 
the signal beam. Therefore, within an angle defined by the separation of 
diffraction orders, k cx - -k, x , k ey - -k gy , and 

E c (x, y, z) = F jj dk 8x dk 8y 

x exp (-ik« • r) |Y(k 8X ,k 8y ) | 2 (k ax ,k 8y ) (174) 

where the propagation vector, k,', is slightly shifted relative to k s caused 
by the frequency difference, 20, between the signal and generated beams. E c 
is an essentially counterpropagating, high-quality conjugate of the input 
field pattern if 7 is a constant when is nonzero. To a good approximation, 
this condition will be met if the spread of the incident wave vectors is small 
compared to the acceptance angle of the waveguide. A constant F requires 
identical waveguides and pump beams that look like counterpropagating plane 
waves, not necessarily at normal incidence. The transverse polarization of 
each waveguide mode in the array and all input beams must be aligned to 
prevent excitation of orthogonal polarization components. 

To quantify the conditions for a high-quality conjugate image, consider, as an 
example, an array where each waveguide element has a Gaussian profile with a 
beam waist of w 0 - lOjrk' 1 . Using Equations 156 and 174 limits the spread of k 
vectors to ~10 mrad for beams no more than 175 mrad (10 deg) from normal 
incidence. The angular separation of array diffraction orders is increased by 
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decreasing the interelement spacing (Ref. 52). A 3« 0 spacing yields an 
angular separation of -70 mrad. The required number of elements in the array 
depends on the length scale of phase-front and amplitude variations in the 
signal beam. The array must be large enough to generate a diffraction-limited 
spot size equal to the inverse of twice the highest spatial frequency in the 
signal beam (Ref. 53). In this example, 1-mrad resolution requires ~2 x 10 4 
elements. For l-/im-wave length light, these numbers correspond to 5-/im 
waveguides on centers of ^50 /im, yielding a total array area of <0.5 cm 2 . A 
signal beam generated by an object of 1 mm 2 in area, located 10 cm from the 
array and containing =100 emitting area elements, would be imaged by the 
backpropagating conjugate beam. This image is also replicated in the multiple 
orders of the array output that fall within the acceptance angle, -AjtOcwq)' 1 , 
of the waveguides. 

If the pump waves are plane waves but are not counterpropagating, then an 
additional phase between waveguides is introduced. For instance, if pump 
beam 1 is not at normal incidence and pump beam 2 is, then the amplitude 
coefficient A^ will have the phase term exp[i(mk lx x 0 + nk ly y 0 ] , as can be seen 
by examining Equation 166. By absorbing the phases into the amplitude 
coefficients, 


^cmn ^cmn ®^P 


- (mk lx x 0 + nk ly y 0 ) 


(175) 


•^amn ^amn ®^P 


i <mk lx x 0 + nk ly y 0 ) 


(176) 


the coupled equations describing the 4-wm interaction in a single-mode 
waveguide return to their original form, but the phase relation between the 
signal and generated beams is changed. This changes the location of the 6 
function that occurs when the summation over waveguides is made in Equa¬ 
tion 173. Now k cll - -k gx + k lx and k cy - -k sy + k ly , changing the propagation 
direction of the generated field but not the fidelity of the replication. The 
conjugate beam generated in a multimode waveguide by 4-wm remains counter- 
propagating when a pump beam angle of incidence is changed (Ref. 11). 
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In conclusion, it has been shown that nearly degenerate 4-wm in an array of 
identical waveguides can yield a high-quality conjugate bean. However, the 
range of wave vectors present in the signal bean must be snail enough so that 
there is nearly equal coupling into the waveguides. Stringent phase-matching 
conditions that apply to the 4-wm interaction in bulk media are lifted by the 
ability to couple light into a waveguide mode over an acceptance angle. Beam 
steering of the generated beam can be achieved by changing the propagation 
direction of one of the pump beams, without destroying the 4-wm interaction 
and while maintaining the spatial coherence of the generated field pattern. 
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6.0 EXPERIMENTS ON OPTICAL FIBERS AND LASER DIODES 


Experimentally testing much of the theoretical framework developed in Sections 
3.0 to 5.0 is impractical given the current minimal availability of complex 
nonlinear waveguiding structures. In this program, two key points have been 
illustrated using a nonpolarization-preserving, single-mode optical fiber and 
GaAs-AiGAs heterostructure and quantum well laser diodes. The first point is 
the importance of removing degeneracy between the waveguide propagating modes 
(Ref. 11). A nonpolarization preserving single-mode optical fiber actually 
contains two degenerate, orthogonal, linearly polarized modes in the weak- 
guiding limit (Ref. 29). The experiments show that polarization scrambling 
between the two modes can be as strong a nonlinear effect as phase modulation. 
The second set of experiments, using the laser diodes, measures the amplitudes 
and relative phases of the output beams generated in the nearly degenerate 
4-wm interaction. The experiments show that, in the weak probe limit, the 
nonlinear interaction in a laser diode acts as a phase modulation of the laser 
diode and that the generated conjugate beam acts to minimize the free carrier 
population pulsations responsible for the third-order nonlinearity. 


EBEIH til*!*•* f3K3QHiH3 il 
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Nonpolarization-preserving, single-mode optical fibers are an example of a 
weakly guiding nonlinear waveguide with degenerate modes. The circularly 
symmetric core leads to two orthogonally polarized modes with identical propa¬ 
gation constants. Small inhomogeneities or temperature and stress effects can 
cause coupling of the two modes even in the linear regime. The purpose of 
this experiment is to show that the fiber nonlinearity causes further modifi¬ 
cation to the polarization of a beam after propagating through the fiber. A 
schematic of the experimental apparatus is shown in Figure 5. A continuous- 
wave (cw) mode-locked Nd:YLF laser, producing =72-ps pulses, is used. The 
beam is passed through a quarter-wave plate and two polarizing beamsplitters 
to produce temporally-spaced, orthogonally polarized pulses which could be 
injected into a 215-m-long silica fiber. A variable beamsplitter is used to 
adjust input intensity. The light transmitted through the fiber was passed 
through a polarizer before detection with a Si-photodiode. 
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Figure 5. Schematic of the experimental configuration Investigating polarization effects In single-mode fibers. 








Initially, the s-polarized path through the beamsplitters was blocked and the 
p-polarized input was detected as a function of power. The polarizer was 
rotated to determine the intensity ratio at the orthogonal angles of maximum 
and minimum transmitted intensity. Figure 6 shows the behavior for two 
consecutive runs. At low intensities the ratio is large. An infinite ratio 
signifies linear polarization. As the intensity increases, the polarization 
ratio tends to oscillate. While the general oscillation trend is clear, there 
are clear differences between the two consecutive runs that are plotted. The 
s-polarized light displays a similar intensity dependence, as shown in 
Figure 7. Note that this third displayed p-polarized data run differs 
slightly from the two plotted in Figure 6. 

In addition to the polarization ratio, the angle of the polarization ellipse 
is important. Figure 8 shows the angle of the polarizer at maximum transmis¬ 
sion intensity, the major axis of the polarization ellipse, for the two 
orthogonal input polarizations as a function of laser power. The reference 
plane is laboratory horizontal and the low-power angle is arbitrary due to 
fiber inhomogeneities so that the relative rotation with increasing power is 
the important measured parameter. The oscillating features of the polariza¬ 
tion ratio data are somewhat mimicked in the ellipse rotation data. Note, 
that the orthogonal polarization relation is not maintained at many power 
levels. 

The importance of the scrambling of polarization on the interpretation of 
other measurements is illustrated by a second set of experiments. The 
modified experimental apparatus is shown in Figure 9. The output optical 
train from the fiber now includes a Mach-Zehnder interferometer. One arm of 
the interferometer has a variable delay so that the delayed s-polarized pulse 
in one arm can be made to temporally overlap with the p-polarized pulse of the 
second arm. In this case, three pulses will arrive at the detector, with the 
output displayed using a sampling oscilloscope. Without polarization scram¬ 
bling there would be no interference between the overlapping s and p pulses. 
However, interference effects can be shown with small displacements of the 
optical delay line. Two cases of relatively large interference effects at 
different power levels are shown in Figures 10 and 11. The interference is 
quite strong at 100 mW/pulse less so at 300 mW/pulse. Referring back to 
Figure 7, the polarization ratio is closest to one near 100 and 300 mW and, 
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POLARIZATION RATIO 


P-POLARIZED INPUT 



Figure 6. The polarization ratio as a function of output power from the 
fiber for the p-polarized input light. 
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Figure 7. The polarization ratio, the maximum to minimum transmitted 

intensity ratio determined by rotating a polarizer, as a function of 
input power. 
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ELLIPSE ROTATION (deg) 


S-POLARIZED INPUT 



TOTAL FIBER OUTPUT (mW) 


Figure 8. The ellipse rotation as a function of output power for 
orthogonally polarized inputs. 








scrambling on interference measurements. 
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Figure 10. Oscilloscope trace of the fast photodiode output using the 
apparatus shown in Figure 9. Average power through the 
fiber is 100 mW. 
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Figure 11. Oscilloscope trace of the fast photodiode output using the 
apparatus shown in Figure 9. Average power through the 
fiber is 300 mW. 
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using Figure 8, the ellipse rotation angle is more nearly orthogonal at the 
higher power level. This explains the relatively strong interference effects 
at both power levels and the greater interference at 100 mW/pulse. 

Polarization changes obscure phase changes due to the degeneracy of the two 
orthogonal modes. Clearly, such nonlinear polarization coupling of modes 
would destroy the fidelity of a PC beam. Deciphering these effects is further 
hindered by the sensitivity of the polarization behavior to external factors 
such as temperature and stress. These types of effects account for the 
variability between the data of Figures 6 and 7. Though displaying 
interesting behavior, interpretation of effects in basic optical fibers is 
less straightforward than the highly nondegenerate slab geometry of 
semiconductor waveguides. 

6.3 LASE R .DI P B£ MMIS 

Four-wave mixing experiments have been performed in single-mode laser diodes 
(Refs. 17, 54 and 55), traveling wave amplifiers (Refs. 40 and 56), and broad- 
area Fabry-Perot devices (Ref. 43). The traveling vave amplifier interaction 
has also been theoretically modeled (Refs. 16 and 41). To date, the 
measurements of 4-wm have concentrated on the amplitude of the amplified 
input, or probe, beam and generated 4-wm signal with little investigation of 
the phase relations between the interacting optical beams. In this work, 
nearly-degenerate 4-wm measurements have been performed in the weak injection 
limit of a laser diode. In this regime, the optical phase of the counter- 
propagating pump beams (^ 0 ). the amplified input signal (<p x ), and conjugate 
4-wm signal ($> 2 ) , satisfy the relation, 


<p 2 ± n = 2<p 0 - <pj 


(177) 


Further, the amplitudes of the electric fields E : and E 2 , corresponding to the 
amplified input and generated 4-wm mixing fields, respectively, are nearly 
equal. This is equivalent to saying that injecting a near-resonant optical 
signal, offset by frequency Q, induces a phase modulation of the laser diode 
at 0. Because of nonlinear interaction results from modulations of the 
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carrier density (Ref. 41), this phase relation means that the generated 4-wm 
beam acts to minimize the density modulations induced by the amplified, 
injected probe. 

A schematic of the experimental apparatus is shown in Figure 12. One laser 
diode acts as a probe. It is optically isolated against back-reflected light 
from any components in the optical train. The output is then passed through 
an acousto-optic modulator where a small fraction, on the order of 10 percent, 
of the beam is deflected and shifted in frequency by 80 MHz. The main beam, 
at an optical frequency i/ x , passes through a variable attenuator and a second 
optical isolator to the measurement apparatus. The shifted beam is injected 
into a second laser diode, operating at a frequency v 0 , slightly shifted from 
j/j by an amount O. The output from the second laser diode is partially 
reflected off of the variable beamsplitter so as to become collinear with the 
v x beam. Both laser diodes are temperature and current stabilized. To bring 
them to nearly degenerate operating frequencies, the temperature is adjusted 
for gross frequency changes and then the current is adjusted for fine tuning. 
Three different laser diodes were used in the measurements. Two are commerci¬ 
ally available GaAs-ASGaAs heterostructure lasers, the Sharp LT015MD (with an 
output power of 30 mW and operating wavelength near 830 nm), and the Hitachi 
HL 7838G (operated at =10 mW and 784 nm). The third type of laser was a GaAs- 
AiGaAs quantum well laser operating near 830 nm and 5-mW output power. This 
laser had a Fabry-Perot cavity length of 600 pm, more than twice as long as 
the commercially available diodes, which gave it greater frequency stability 
(Ref. 57). 

To measure the spectrum of the emitted light, two spectrum analyzers were 
used. The first was a Newport SR-240 high-finesse scanning Fabry-Perot with 
detector. This instrument had a free spectral range of 8 THz and finesse of 
20,000, giving an optical frequency resolution of 400 MHz. Alternatively, the 
laser diode output was detected by a fast photodiode, 3-dB point at 2 GHz, and 
the photodiode signal was displayed on a radio frequency (RF) spectrum 
analyzer. Because the photodiode is a square-law detector, its output will 
have a spectrum corresponding to the difference in frequencies between various 
optical frequency components. 
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LD = LASER DIODE 


Q = LENS 

—ft- = FARADAY OPTICAL ISOLATOR 
A.O. = ACOUSTO-OPTIC MODULATOR (80 MHz) 
| = MIRROR 
J = BEAMSPLITTER 
=. PHOTODIODE 

O.S.A. = OPTICAL SPECTRUM ANALYZER 
RF S.A. = RF SPECTRUM ANALYZER 


Figure 12. Schematic of the experimental arrangement for 4-wm in laser 
diodes. 

Detection of four optical frequency components can be expected. They are the 
two diode operating frequencies, u 0 and i/ lt the acousto-optically shifted 
frequency, v x + 80 MHz, and the 4-wm frequency, 2v 0 - (•'i + 80 MHz). However, 
80 MHz is an insufficient frequency shift to be resolved by the Fabry-Perot. 
With the RF spectrum analyzer, there will be potentially many mixing 
components, but assuming that the optical frequency components at t/ 0 and v x 
are strong compared to the others, the primary frequency components should be 
at 80 MHz, O, 2 0 ± 80 MHz, and n ± 80 MHz. The plus or minus sign will 
depend on the magnitude of with respect to i/ 0 . The first three radio 
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frequencies result from the mixing of the three output frequencies from laser 
diode #2 with u 1 , while the fourth is the beat frequency between i/ 0 and its 
two sideband frequencies. 

Figure 13 displays the output of the optical spectrum analyzer when the two 
laser diodes are tuned to near degeneracy. The main peak is at u 0 and the 
larger side-peak is primarily at i/ x . Oppositely shifted from v Q is the 
generated nearly-degenerate 4-wm signal. If the t/ x beam is blocked just 
before the variable attenuator, only the output from laser diode #2 is 
displayed. This is shown in Figure 14. Now the amplified version of the beam 
injected at v x + 80 MHz is shown to the left of the v Q signal. Note that it 
is at approximately the same magnitude as the 4-wm signal. 

When the fast photodiode and RF spectrum analyzer are used, a typical RF 
spectrum is shown in Figure 15. The frequency range that could be measured 
was limited by a 500-MHz bandwidth amplifier that followed the photodiode. 
There are three RF features that result from the mixing of optical frequen¬ 
cies. The two small peaks near zero frequency are noise artifacts. An 
instrument resolution limited feature at 80 MHz is generated by v x and the 
amplified probe signal at v x + 80 MHz. This signal does not result from 
reflections off of optical components. This was verified by tuning the 
current of laser diode #2 which moved the operating frequencies. Unless u x 
and v 0 were nearly degenerate, the 80-MHz signal disappeared. The tuning 
range will be discussed further. The second, and largest, feature in the RF 
spectrum is due to the mixing of v 0 and i/j. It is broadened due to the 
inherent frequency jitter of the free-running diodes. Finally, at higher 
frequency, there is a second broadened feature at a frequency corresponding to 
2 0 to 80 MHz. Therefore, it is the mixing signal between i/j and the 
generated 4-wm frequency. 

Significantly, there is no RF feature at the difference frequency fi - 80 MHz, 
though there are certainly optical frequency components with this spacing. 

This is even more surprising given that the strength of the optical power 
detected at v 0 is more than an order of magnitude larger than at . Refer¬ 
ring back to Figure 14, note again that the two sideband components are of 
nearly equal amplitude. The photodiode signal, measuring optical power rather 
than field amplitude, shows no power modulation at fi - 80 MHz. This point 
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OPTICAL POWER (ARB. UNITS) 



Figure 13. Optical spectrum analyzer trace using the apparatus shown in 
Figure 12. The upper trace is an offset, lOOx expansion of the 
lower trace. 
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Figure 14. Optical spectrum analyzer trace using the apparatus shown in 

Figure 12, except that the direct optical path from laser diode #1 
to the beamsplitter has been blocked. The upper trace is an 
offset, lOOx expansion of the lower trace. 
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Figure 15. Radio frequency spectrum analyzer trace using the apparatus shown 
in Figure 12. 

requires that the optical sidebands represent a phase modulation of v 0 and 
that the phases of the three frequency components generated by laser diode #2 
be related according to Equation 177. The fact that the RF feature at 80 MHz 
is resolution limited means that the amplified input signal remains injection- 
locked, while the broadening of the other two features shows that the two 
lasers remain unlocked. 

It is now generally accepted that the 4-wm signal in laser diodes results from 
free-carrier population pulsations in single-mode devices and gratings in 
broad-area devices (Refs. 41 and 43). The pulsations and gratings are assumed 
to be set up by modulation of the carrier density at 0 due to the injected 
beam. The phase relation between the three optical beams means that the 
nonlinear interaction acts to limit the population pulsations induced by the 
injected signal. Using Agrawal's theory (Ref. 41), if the amplified injected 
signal were absolutely equal in amplitude, there would be no free-carrier or 
refractive index modulations. The generated beam acts to cancel the intensity 
modulations experienced in the laser diode and keep the carrier pulsations to 
a minimum. Previous work has shown that the generated 4-wm signal saturates 
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with increasing injected signal power (Ref. 54). This implies that the 
amplified injected signal becomes less efficient at modulating the carrier 
density as it tends toward saturation. 

A requirement of an optically-induced real population pulsation or grating 
through 4-wm is that the beams involved in the mixing have collinear polariza¬ 
tion components. There is no field mixing with orthogonal polarizations. 

Using a half-wave plate inserted in the i/ x - 80 MHz beam path, the amplified 
injected beam signal and 4-wm signal were verified and resulted only from the 
polarization component injected parallel to the circulating laser field 
polarization in laser diode #2. As the half-wave plate was rotated, the 
amplified injected signal and 4-wm signal scale proportional to the intensity 
of the component parallel to the circulating laser-field polarization. The 
orthogonal (vertical) polarization was neither amplified nor did it induce a 
signal in the parallel (horizontal) polarization direction. 

There are several interesting effects which appear as the frequency of the 
signal injected into the diode is varied. Figures 16-18, along with Figure 15 
are a sequence of RF spectra shoving changes as the injection laser is scanned 
through the oscillation frequency of the second laser diode. In all cases, 
the injected signal is very weak so that the 4-wm signal is «10' 3 of the main 
signal strength. Continuing the convention, the laser diode #1 is at v lt the 
acousto-optic modulator generates a beam at + 80 MHz and the second laser 
diode operates at i/ 0 . Figure 15 is the spectrum with j/j tuned well below i/ 0 . 

As the frequency of is increased, the v 0 - v lt and 2(v 0 - Vj) - 80 MHz 
features shift to lower frequencies until the injected signal locks the second 
laser diode. An injection-locked spectrum is shown in Figure 16. There is a 
single, sharp feature at 80 MHz, the acousto-optic frequency shift. The 
feature sits atop a pedestal that is broadened by the phase and amplitude 
noise of the two diodes. Continuing to increase the frequency of i/ 1( now >i/ 0 , 
the injection-locking range is less than «50 MHz with this weak signal, so 
that when t/ x - v 0 there are simultaneously narrow and broad features, at 
80 MHz, as shown in Figure 17. The narrow feature is the amplified-injected 
signal and the broad feature is the symmetrically offset 4-wm signal with the 
frequency inverted. There are two weak features at 160 and 240 MHz which are 
electronic noise generated by the acousto-optic modulator power supply. 
Finally, as vj is increased further, the - u 0 peak and 4-wm peaks move to 
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Figure 16. Radio frequency spectrum analyzer trace using the apparatus shown 
in Figure 12 when the two laser diodes are phase locked. 
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Figure 17. Radio frequency spectrum analyzer trace using the apparatus shown 
in Figure 12 when the two laser diodes have a frequency difference 
of -80 MHz. 
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Figure 18. Radio frequency spectrum analyzer trace using the apparatus shown 
in Figure 12 when laser diode #2 operates at a higher frequency 
than laser diode #1. 

higher radio frequencies, again reflecting the RF spectra display of the 
absolute frequency difference. In this case, however, the 4-wm signal is now 
at 2(»' 1 - u 0 ) + 80 MHz, as shown in Figure 18. Except over the narrow 
injection-locked frequency range, the spectra are qualitatively similar. 

A second series of spectra, taken with the optical spectrum analyzer, are 
shown in Figures 19-23. The injected power is approximately a factor of 5 
larger than in Figure 15 in this series. Again, frequency increases from left 
to right and only the output from laser diode #2 reaches the detector. In 
Figure 19, the injected signal is the small peak «8 GHz below the main 
feature. Also, there is an even smaller, oppositely-shifted 4-wm signal. 

There are also symmetric shoulders on the main peak that correspond to the 
relaxation oscillation resonance (Rsfs. 17 and >8). The resonance results 
from the photon-free carrier coupling and is at a frequency between the 
inverse laser cavity photon lifetime, on the order of a few picoseconds, and 
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Figure 19. Optical spectrum analyzer trace of the output from laser diode #2 
when the injection frequency is offset by ~8 GHz. The upper trace 
is a 50x expansion of the lower trace. 
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Figure 20. Optical spectrum analyzer trace of the output from laser diode #2 
when the injection frequency is offset by the relaxation resonance 
frequency. The upper trace is a 10x expansion of the lower trace. 
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Figure 22. Optical spectrum analyzer trace of the output from laser diode #2 
when the injection frequency is offset by ~1 GHz. The upper trace 
is a lOx expansion of the lower trace. 
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Figure 23. Optical spectrum analyzer trace of the output from laser diode #2 
when it is phase locked. The upper trace is a 10x expansion of the 
lower trace. 


inverse free carrier lifetime, on the order of a nanosecond. When the injec¬ 
tion frequency is offset from the free-running frequency by greater than the 
relaxation frequency, as shown in Figure 19, the 4-wm interaction and the 
amplification of the injected signal are quite weak. At the relaxation 
resonance, these signals are nearly a factor of 50 larger, as shown in 
Figure 20. Peaks at twice the relaxation frequency, due to the relaxation 
peaks of the injected signal mixing with the sidebands of laser diode #2, are 
also visible. The signals are reduced by somewhat more than a factor of 3 at 
half the relaxation frequency (Fig. 21) and only become as strong again when 
the offset is reduced to =1 GHz, as shown in Figure 22. Also note that there 
is some asymmetry in the amplitude of the two sidebands, indicating that 
amplitude as well as phase modulation is playing a role in the laser diode 
operation. Figure 23 shows injection-locked operation. These results are 
consistent with those previously reported (Ref. 17). 
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A final series of spectra taken with the optical spectrum analyzer are shown 
in Figures 24-28. These were made using the Sharp laser diodes, which had a 
much stronger relaxation oscillation than the Hitachi diodes measured in 
Figures 19-23. For these spectra, the injected signal was quite strong, 
approximately an order of magnitude larger than the previous series. The 
free-running output of laser diode y/2 is shown in Figure 24, and a weakly 
amplified input signal well offset is shown in Figure 25. 

As i/ x is tuned toward the relaxation resonance peak, the frequency spectra 
change dramatically. Representative spectra are shown in Figures 26-28. The 
offsets are = -6, -3, and +3 GHz, respectively. Figures 27 and 28 show strong 
evidence of a set of closely-spaced or continuous-frequency components over a 
range of the spectra. This could be evidence of coherence collapse and/or 
chaotic phenomena, but verification requires a more detailed study with higher 
frequency resolution than was available in the current program. 

The spectra identify several key features. First, there are at least three 
operating ranges of injection intensity. At low levels, the injected signal 
acts as a phase modulator generating symmetric sidebands. At intermediate 
levels, the sidebands become asymmetric while at high levels spectral broaden¬ 
ing and multiple asymmetric features dominate the optical spectrum. A key 
parameter in determining signal strength is the relaxation resonance. This is 
to be expected given the basis of the 4-wm signal as the coupling between the 
photon and free carrier densities. Finally, it should be noted the two lasers 
could never be stably injection-locked by the 4-wm sidebands. When the two 
lasers showed locking behavior, the output frequencies of the two laser diodes 
were always degenerate. 
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Figure 24. Optical spectrum analyzer trace showing the free-running output of 
the Sharp LT015 laser diode with its relatively strong resonance 
sidebands at ~ ±3 GHz. 
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Figure 25. Optical spectrum analyzer trace showing the LT015 laser diode 
under strong optical injection ~10 GHz below the main peak. 
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Figure 26. Optical spectrum analyzer trace showing the LT015 laser diode 
under strong optical injection ~6 GHz below the free-running 
frequency. 
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Figure 28. Optical spectrum analyzer trace showing the LT015 laser diode 
under strong optical injection ~3 GHz above the free-running 
frequency. 
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7.0 CONCLUSIONS 


This program has investigated 4-wm in nonlinear waveguides, with a 
concentration on generation of a conjugate output beam in materials that 
exhibit both a strong nonlinearity and resonance with linear gain. The 
coupled-mode theory was used to develop a general formalism to describe 
nonlinear optical interactions in a waveguide. Under a wide range of 
experimentally and technologically relevant configurations, the basic result 
can be summarized: 

1 y, 

= ± 1 2 ^ <l> exp(iP a z) // dxdy P 111 * (r,o>) • i|; a (x,y) (178) 

Here, the ± refers to propagation along the ±z (waveguide) direction, the 
l(r a (x,y) and /J a are the linear mode transverse field patterns and propagation 
constants, respectively, and the A a are the amplitude coefficients for each 
mode, which vary due to the nonlinearity. The modes are used to describe the 
electric field patterns: 

E(r,t) = £ A a (z)ijr a (x,y) exp [i (|J a z-tot) ] (179) 

a 

The mode amplitudes become coupled through the effect of the nonlinear 
polarization: 

^(r, &> 0 ) = x i2) (*/<•> = «1 + <*> 2 ) : E(r,u 1 ) E(r,o> 2 ) 

+ X <3> (r,o> = + u >2 + u 3 ) : E(r,u 1 ) E(r,o 2 ) E(r,o> 3 ) 

+ (180) 

An important feature of Equation 178 is its similarity to the coupled-wave 
formalism which has been used to successfully model nonlinear optical inter¬ 
actions in bulk media. Nonlinear phenomena in waveguides, such as 4-wm, can 
be modeled similarly to the bulk interaction, with the overlap integral of 
Equation 178 accounting for the spatial confinement effects. An important 
requirement for the use of a single amplitude coefficient, A a , with each 
guided mode is that all fields involved in the nonlinear interaction can be 
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described by Che propagating modes of a weakly guiding system. Though the 
coupled-mode theory can handle situations where all modes are not confined and 
the weakly guiding approximation does not hold, these situations will almost 
always tend toward reduced fidelity in replication of an input beam through 
the 4-wm process. This result, in combination with results of Hellwarth's 
earlier analysis (Ref. 11), puts strict limits on the speed and image fidelity 
of PC generation by 4-wm in nonlinear waveguides. 

Of particular interest in this study are materials where the nonlinearity 
results from density gradients of mobile particles, such as semiconductors 
under current injection. In these materials, carrier diffusion acts to limit 
gradients of high-spatial frequency. Semiconductor waveguides generally have 
rectangular cross sections and transverse mode patterns which are, in general, 
described by two transverse spatial frequencies (Ref. 50). Nonlinear optical 
coupling will only occur between those modes where the spatial frequency 
difference is, approximately, less than the inverse diffusion length. This 
limits the divergence angle of an optical beam, or the angle between two 
incident beams, to a few degrees in a typical III-V semiconductor material. 

The experimental results of Lucente, et al. (Ref. 43) are consistent with this 
analysis. 

To overcome these limitations in multimode waveguides, the use of an array of 
uncoupled, identical, single-mode waveguides was proposed. In addition to 
offering the potential of high-fidelity phase conjugation, the array can 
produce conjugate wavefronts whose propagation direction can be controlled by 
the pump beams. Unfortunately, a suitable array cannot be obtained to test 
the predictions of the model. In fact, understanding of 4-wm using two or 
three single-mode waveguides has yet to be achieved. 

A particularly important waveguide is the laser diode. A variety of nonlinear 
phenomena, from saturation to chaotic behavior, have been observed. Specific¬ 
ally, 4-wm in laser diodes and traveling wave amplifiers has been investigated 
experimentally (Refs. 17, 43, 54 and 55) and theoretically (Refs. 16 and 41). 
Past work has shown large nonlinear amplification due to optically-generated, 
free-carrier density pulsations or gratings. The nearly degenerate 4-wm 
interaction has been shown to cause the laser diode to respond as a phase 
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modulated device, in the limit of a weak injection signal. In effect, the 
laser diode acts to minimize carrier modulations by generating a new frequency 
components. 

This result, when compared with previously published results on 4-wm in laser 
diodes, raises several tantalizing questions. First, all previous results 
were made outside of the weak injection signal limit and show asymmetric 
sidebands on an optical spectrum analyzer. Because no RF spectra have been 
taken, carrier response in this regime has not been determined. The influence 
of two or n input beams on the diode output has not been measured, nor has the 
influence of the transverse coupling that exists, for example, in a diode 
array. Further, the nonlinear interaction medium is a Fabry-Perot. Placing a 
material with a strong nonlinear optical response within an optical cavity 
couples the amplitudes and phases of the circulating optical fields. While 
the importance of the optical feedback on the 4-wm has been recognized, the 
specific effects, particularly on phase relations between interacting beams, 
are not clearly understood. Recently, it has been demonstrated that a diode 
laser can be mode-locked when it is coupled to a linear, empty external cavity 
which provides direct feedback to the laser.* This implies strong nonlinear 
interactions between the laser modes under self-feedback. Because of the 
importance of optical phase relations to problems of technological interest, 
including the generation of a coherent optical output by arrays of laser 
diodes, further investigation of the basic physics of multiwave mixing in 
active nonlinear waveguides is an important task. 


* G. Yao, P. Wang, Y.C. Chen, C.M. Harding, R.J. Dalby, R.G. Waters, J.M. Liu, 
and K.K. Lee, "Self-Starting Additive-Pulse Mode Locking in Semiconductor 
Lasers with Linear External Cavity," submitted to Opt. Lett. 
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APPENDIX A 

DOT PRODUCT OTHOGONALITY RELATIONS 


A simple set of dot-product orthogonality relations can be derived for CASES 2 
and 3, when the propagation axis is an axis of symmetry for e. This follows 
from the identity 


V • (E a x H b ) 



+ 


i(i> 

c 




(A-l) 


Integrating the equation over the transverse cross-sectional area and using 
Equations 52 and 53 in Subsection 3.3.2, the left-hand side of the equation 
vanishes for * fa. This leaves 


f J dxdy (i|f B * e • V b * <p, 


• <p b ) = 6 


o> 


a,b 


(A-2) 


Using the relations between the forward- and backward-propagating modes in 
Subsection 3.3.1 yields 


JJdxdy (i|r aT • € t • iJr bT + <p { 


r«Pbx) * (*..b * 


(A-3) 


j‘]‘dxdy(e z2 tjf az i|r b2 + <p aT • <p bT ) = (6 a , b - 6 a . b ) 


(A-4) 


Equations A-3 and A-4 hold for both real and complex e(w) so long as z is a 
symmetry axis for e. No such simple relations are possible when z is not a 
symmetry axis. 

It is important to note from Equations A-3 and A-4 that, in general, 


JJdxdyijr a • c • ^r b * 0 and J* J dxdy\|r aT • e T • <|r bT * 0 (A-5) 

JJdxdy<p a • <p b * 0 and Jjdxdy<p aT • <p w * 0 ( A * 6 ) 
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When the transverse variation of the mode fields is slow on the scale of & 1 , 
then the z-component will be small relative to the transverse part of the 
fields. In this case, the integrated transverse dot product of different 
modes will be approximately zero. Alternatively, this term will be zero in 
specific geometries even though the z-component exists, such as the slab 
geometry where there is no dielectric variation along a transverse direction. 
For systems with cylindrical symmetry in general, the relations in Equations 
A-3 and A-4 complement the cross-product relations derived in Subsection 




APPENDIX B 


ITERATIVE PROCEDURE FOR E WHEN P 1 ^ 4 * 0 

When P 1 ^ 4 ¥■ 0, Equations 62 to 65 have to be used to expand the nonlinear 

field components in a waveguide. The mode coefficients, A a (z), cannot be 
solved directly through the coupled-mode equations. Instead, they can be 
solved self-consistently through an iterative procedure outlined below. 

1. Initially set P 1 ^ - 0. 

2. Write 


Er = E A a (z) t|r aT (x,y)exp(i(J a z) (B-l) 

a 

E z = E A a (z)i|r az (x,y)exp(ip a 2 r) - pf* (B-2) 

a € zz 

In this step, A a (z) are unknown coefficients to be solved later, but an 
explicitly known P 1 ^ 4 given by the initial condition or calculated from the 
last iteration is used. 

NL NL 

3. Using Equation 15, express P T and P z for the current iteration in 
terms of A a (z). 

4. Write coupled-mode equations using Equation 81 and solve for A a (z). 

5. Calculate P 1 ? 4 and Pf explicitly using expressions obtained in 
Step 3. Calculate and E z using the newly calculated values of 
A a (z) and P 1 ? 4 . 
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Compare the new values of E and P 1 * 1 * with those of the last 
iteration. 

a. If significant difference exists, go to Step 2 using the newly 
calculated value of for the input. 

b. Otherwise, a self-consistent solution for E and P*™* is obtained. 






